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ABSTRACT 

A physical description of the fundamental plane of Galactic globular clusters is 
developed which explains all empirical trends and correlations in a large number 
of cluster observables and provides a small but complete set of truly independent 
constraints on theories of cluster formation and evolution in the Milky Way. 

In a very good first approximation, the internal structures and dynamics of 
Galactic globulars are described by single-mass, isotropic King (1966) models and 
are thus characterized by four nominally independent physical parameters. Within 
this theoretical framework, it is shown that (1) 39 regular (non-core-collapsed) 
globulars with measured central velocity dispersions all share a common core 
mass-to-light ratio, Ty,o = 1-45 M@ Lq , and (2) 109 regular globulars both with 
and without direct observations of Tv,o show a very strong correlation between 
global binding energy and total luminosity, regulated by Galactocentric position: 
Eb = 7.2 x 10 39 erg (L/LQ) 2 - 05 (r gc /8kpc) -a4 . The observational scatter about either 
of these two constraints can be attributed fully to random measurement errors, making 
them the defining equations of a fundamental plane (FP) to which real clusters are 
confined in the larger, four-dimensional parameter space of general King models. 
They are shown to underlie a pair of bivariate correlations first found, and used 
to argue for the existence of a globular cluster FP, by Djorgovski (1995). A third, 
weaker correlation, between clusters' total luminosities and King-model concentration 
parameters, is related to the (non-random) distribution of globulars on the FP. With 
L, Tv,o, Eb, and the central concentration c thus chosen as the four physical quantities 
that define any single globular cluster, the FP equations for T^o and Eb(L,r gc ) are 
used to derive expressions for any other observable in terms of only L, r gc , and c. 
Results are obtained for generic King models and applied specifically to the globular 
cluster system of the Milky Way. 
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1. Introduction 

The properties of globular clusters in the Milky Way offer empirical constraints not only 
on their own formation and evolution, but also on the history of our Galaxy as a whole and, to 
a potentially large extent, on the formation of generic star clusters — a key element of the star 
formation process itself. The many observed correlations between various internal structural and 
dynamical parameters of the clusters (core and half-light radii, surface brightnesses, velocity 
dispersions, etc.) that have been found by several authors (e.g., Brosche & Lentes 1984; Kormendy 
1985; van den Bergh, Morbey, & Pazder 1991; van den Bergh 1994; Djorgovski & Meylan 1994; 
Djorgovski 1995; Bellazzini et al. 1996; Murray Sz Lin 1992) constitute an important set of such 
constraints. 

The apparently large number of these correlations (see especially Djorgovski & Meylan 1994) 
can be somewhat misleading, however, in that it is impossible for all of them to be independent. 
This follows from the fact that Galactic globulars are structurally and dynamically exceedingly 
simple, being very well described by single-mass, isotropic King (1966) models. Since these 
models are defined fully by the specification of just four independent physical parameters, it is 
not unreasonable to expect that the same should be true, at least in a first-order approximation, 
of real clusters; but there can then be at most three truly distinct correlations between their 
various intrinsic properties (although these may be supplemented, in general, by dependences on 
"external" factors such as Galactocentric position or cluster metallicity) . 

It is clearly of interest, then, to systematically reduce the globular cluster data in the Milky 
Way to the minimum number of independent physical relationships required for a complete 
description of the observed properties of the ensemble. The establishment of such a compact set 
of empirical constraints — which should serve better as a target for theories of cluster formation 
and evolution — is the main goal of this paper. To reach it, two important results will be used as 
starting points. 

First, the two best globular cluster correlations known to date (i.e., the only two with scatter 
that can be fully accounted for by observational uncertainties) are the bivariate dependences 
presented by Djorgovski (1995) for a subset (roughly one-third) of the Galactic globulars. One 
of these correlations involves only cluster core parameters, while the other includes properties 
measured at their half-light radii. Djorgovski found these correlations by a statistical principal- 
components analysis, and he offered no physical interpretation of them; but he pointed out that 
their very existence suggests that globular clusters may in practice be confined — much like elliptical 
galaxies and other large stellar systems are — to a "fundamental plane" in the higher-dimensional 
space of physical parameters available to them in principle. 

Second, an almost completely unexamined (but, as it happens, very closely related) issue has 
to do with the dependence of binding energy, E^, on mass (or luminosity) among globular clusters. 
This was studied some time ago by Saito (1979), who also investigated the question for dwarf and 
giant elliptical galaxies. He found that Ej, oc M L5 or so — apparently consistent with very simple 
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virial-theorem arguments — for globulars and giant ellipticals alike. However, his observational 
sample included only ten relatively bright globular clusters. The data required for a calculation of 
binding energy are now available for many more clusters than this and they are obviously of higher 
quality than those available to Saito. A total energy Ef, is the only obvious physical property that 
is not routinely estimated from observations of Galactic globulars; but it is arguably important — if 
for no other reason than to complete a systemization of the data — to clarify the empirical behavior 
of this quantity. 

In what follows, catalogued data are used to estimate binding energies within the context of 
single-mass, isotropic King (1966) models for 109 "regular" (non-core-collapsed) globular clusters 
in the Milky Way. The data and necessary aspects of the models are summarized in the next 
Section and in a pair of Appendices. Section 3 presents the binding energies themselves, and 
also derives core mass-to-light ratios, Tyo> f° r 39 of these clusters with measured central velocity 
dispersions. It is shown that (1) Yy,o is a constant for all clusters, and (2) Ej, correlates tightly 
with total cluster luminosity (or mass), scaling as Ej, oc L 2 05±ai . This is significantly (if not 
surprisingly) different from Saito's (1979) result, and perhaps also from the relation obeyed by 
elliptical galaxies (although this also appears not to have been reviewed in any detail since Saito's 
work) . 

Section 3 also establishes the choice of L, Ej,, Yyo, and the King-model concentration 
parameter c = log (r^/Vo), as the four "independent" physical properties used to define any 
globular cluster. The tight constraints found on Tyfl and Ei{L) then imply that real clusters 
have (at most) only two free parameters, and they do indeed occupy a fundamental plane in the 
theoretical space of King models. The influence of environment on this plane, as reflected in 
cluster Galactocentric radii and metallicities, is also examined in §3; it is shown to be confined to 
a dependence on r gc in the normalization of the E^-L scaling. 

In §4, an orthogonal coordinate system called e-space (analogous to the so-called K-space of 
elliptical galaxies; Bender, Burstein, & Faber 1992) is constructed to obtain a face-on view of the 
fundamental plane. This shows that globulars are not distributed uniformly over it; essentially, the 
parameter c, which describes the internal density profile of a cluster, is also correlated with total 
luminosity or mass (see also van den Bergh 1994; Bellazzini et al. 1996). Alternate representations 
of the fundamental plane, in the form of a number of bivariate correlations between cluster core 
and half-light parameters, are also developed in §4. It is shown that the two correlations of 
Djorgovski (1995) are equivalent to the two constraints Tyo = constant and Ei oc L 2 ' 05 and, thus, 
that Djorgovski identified the same fundamental plane that is developed here. 

Section 5 and Appendix A then use the defining equations of the fundamental plane, along 
with manipulations based on generic properties of King model, to derive a set of simple relations 
that ultimately describe the variation of any globular cluster observable as a function of L, c, 
and r gc only. These may themselves be manipulated to explain any observed correlation between 
any other combination of cluster variables. Thus, the fundamental plane indeed reduces a large 
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array of globular cluster data to a very small set of independent constraints for models of their 
formation and evolution. These are set out in §6, which summarizes the paper. 

It should be noted that the choice of "basic" cluster properties adopted here — total luminosity 
(mass), binding energy, mass-to- light ratio, and concentration parameter — is not a unique 
possibility, nor is it necessary in any rigorous sense. It will be noted, for example, that the scaling 
of binding energy with luminosity is essentially equivalent to other, previously known correlations 
between half-light radii and luminosities (Rh ~ L°; van den Bergh et al. 1991) and between core 
velocity dispersions and luminosities (<To ~ L 0,52 ; e.g., Djorgovski & Meylan 1994). However, 
this does not make a description in terms of E^ any less correct or valid. In fact, this framework 
seems particularly amenable to theoretical studies of cluster evolution and formation both: Many 
simulations of the former follow the evolution of c and Ef, over time, and a discussion of the 
Eb(L,r gc ) dependence in the latter context will be taken up in a subsequent paper (McLaughlin, 
in preparation). 



2. Data and Model Calculations 



It is well known that King's (1966) theoretical sequence of isotropic, lowered isothermal 
spheres of single-mass stars (see also Binney & Tremaine 1987) provides an excellent match to 
the luminosity profiles of most Milky Way globulars and a good description of their internal 
kinematics. This specific model framework is therefore adopted for analyses throughout this 
paper. For self-consistency, and in order to make use of as large an empirical database as possible, 
the complications of velocity anisotropy and multi-mass stellar populations are not considered 
(although theoretical models incorporating these features have, of course, been constructed [Michie 
1963; Gunn & Griffin 1979] and applied to a number of Galactic globulars [see Pryor & Meylan 
1993]). 

Harris (1996) maintains an on-line catalogue^] of King-model parameters and other properties 
of 147 Galactic globular clusters; it is in part a synthesis and updating of earlier compilations by 
Trager, Djorgovski, & King (1993) and Djorgovski (1993a) (see also Trager, King, & Djorgovski 
1995). Among other quantities, the catalogue lists the clusters' heliocentric distances and their 
Galactocentric positions, r gc ; their foreground reddenings, E(B — V) = Ay/3.1, and absolute 
magnitudes, My; metallicities, [Fe/H]; projected half-light (or effective) radii, R^; central V-band 
surface brightnesses, /iv,o (not corrected for extinction in the raw catalogue); scale radii,[] ro, 
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ittp : //phy sun, physics .monaster . ca/Globular .html . The version used in this paper is the revision of 22 June 



The radii ro are referred to by Harris (1996), and by Trager et al. (1993) and Djorgovski (1993), as core radii, 
and they are denoted r c in these catalogues. However, the tabulated values in fact correspond to the King-model 
scale radii ro defined by equation (^), and they are quantitatively different from the projected half-power radii of 
low-concentration clusters. 
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and concentration parameters, c = log (rj/ro) (with rt the tidal radius), derived from fits of King 
(1966) models to cluster surface-brightness profiles; and central luminosity densities, jo (in F-band 
L & pc -3 ), obtained from the measured c, r , and extinction-corrected nv,o (see, e.g., Djorgovski 
1993a, and Appendix A below). 

The King- model structural parameters are not available for six of the globulars in Harris' 
(1996) tables; another cluster (Djorgovski 1) has no measured metallicity; and no central surface 
brightness is given for one more (the Pyxis cluster). These objects were therefore excluded a priori 
from the data set ultimately used here. Of the 139 clusters remaining, 30 are identified (following 
Trager et al. 1993) as having obvious or suspected post-core-collapse (PCC) morphologies, i.e., 
central density cusps. As Trager et al. discuss, such clusters are not particularly well described 
by King models — a concentration parameter of c = 2.5 has generally been assigned to them 
arbitrarily — and their tabulated scale radii ro and core surface brightnesses fiyfl are most likely 
over- and underestimates, respectively. Thus, although most of the calculations described below 
have been performed for the PCC clusters using the data as given, the results should not be taken 
at face value; and while these objects are shown on many of the plots in this paper, they are never 
included in any more detailed analyses (e.g., least-squares fits). 

This leaves 109 regular, King-model clusters as the main focus of attention. Of these, 103 
have Galactocentric radii r gc < 40 kpc; 4 have 65 <^ r gc ^ 100 kpc; and two (the faint Palomar 
4 and AMI) lie at r gc ~ 110-120 kpc. Estimates of the central line-of-sight velocity dispersions, 
<jpfi, are given for roughly one-third (39/109) of this restricted sample, as well as for 17 of the 30 
PCC clusters, by Pryor & Meylan (1993). These additional data are also employed here. 

The recent study of Rosenberg et al. (1998) has been used to replace some of Harris' (1996) 
numbers for the low-luminosity cluster Palomar 12. Specifically, from a King-model fit to star 
counts, Rosenberg et al. derive c = 1.08 and ro = Of 63 (or 3.5 pc, for a heliocentric distance 
of 19.1 kpc), as compared to the c = 1.94 and ro = 0'.20 given by Harris. With nvfl = 20.6, 
according to Harris, this then implies (following the procedure outlined in Appendix A below) a 
central luminosity density jo ~ 32 L Q pc -3 , lower than the previous value by a factor of roughly 
three. After incorporating these few changes for this one cluster, the final database is similar to 
that adopted by Djorgovski & Meylan (1994), and the subsets thereof in Djorgovski (1995) and 
Bellazzini (1998), for their recent investigations of the Milky Way globular cluster system. 

Expositions of observational uncertainties can be found in Trager et al. (1993) and Pryor & 
Meylan (1993), and the situation may be summarized roughly as follows: scale and (projected) 
half-light radii are known to typical precisions of A log ro — ±0.1 dex and A log ~ ±0.1 
dex; central velocity dispersions, to within A log a P) o — ±0.07 dex; central surface brightnesses, 
to A nvfi — ±0.3 mag arcsec -2 ; and concentrations, to Ac ~ ±0.2 dex. Uncertainties in 
any derived quantities follow immediately from these estimates. For example, the central 



luminosity density is obtained from /j,y,o = const. — 2.5 log (jo r o) (eqs. [A2|, JA3||), and thus 



A log jo — ±0.2 dex. Similarly, the average surface brightness within a half-light radius is defined 
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as (pv)h = const. — 2.5 log (L/Rfy (see eq. | A5|), so that A (pv)h — ±0.5 mag arcsec -2 . There are 
also uncertainties in cluster distance moduli — ±0.15 mag may be representative; see the discussion 
at Harris' (1996) website — and thus in Galactocentric distances r gc and total luminosities L (note 
that integrated apparent V magnitudes are generally good to a few hundredths of a magnitude). 
These are generally not dominant in the error budgets discussed here, however, and for the most 
part they can be ignored. 

In order to refer cluster data to isotropic King models, it is necessary first to relate the 
observed velocity dispersion a Pt o to the model scale velocity, a®. This latter quantity is defined 
through the distribution function assumed for the cluster stars, f(E) oc cr^ 3 [exp (E/o-q) — l] , 
and through the dimensionless cluster potential, W(r) = [4>{rt) — <p{r)] /o~q (King 1966; Binney 
& Tremaine 1987). W(r) and a scaled density profile, p/ po vs. r/r , follow from integrating 
Poisson's equation with this distribution function and a specified value of Wo = W(r = 0). With 
rt the radius at which W = p = 0, there is then a one-to-one relation between Wq and the 
concentration parameter c = log (rj/Vo). In the limit of high concentration, the King model is a 
normal isothermal sphere, ctq is the true one-dimensional velocity dispersion inside the cluster, 
and o~ p fi = o"o- For low concentrations c <J 1, however, the outer, non-isothermal parts of a cluster 
influence projected quantities even at its center, and <t Pj o < o~q. This effect is illustrated in the 
upper left-hand panel of Fig. |]. The numerical integrations summarized there have been used 
(given the observed concentrations c) to convert the central dispersions a Pt o to model o"o-values 
for the 39 regular and 17 PCC clusters in Pryor & Meylan (1993). In most cases, this is not a 
large correction: the Pryor & Meylan clusters all have c > 0.75, and 0.8 < o^o/co < 1- It does 
have the effect of adding slightly to the observational uncertainty, however: with Ac = ±0.2 and 
A log <7 Pi o = ±0.07, the typical errorbar on log uq is roughly ±0.09 dex. 

The binding energies of Milky Way globulars are calculated in §3 from the basic definition 
(cf. King 1966; Saito 1979) 
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where M is the total mass and the sign has been chosen to make E\, > for bound clusters. A 
dimensionless binding energy follows directly from this and is a function only of c (through Wo) 
for an isotropic and single-mass King model: 
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The function £ is shown in the upper right panel of Fig. |l|. Over the range 0.5 < c < 2.5 
appropriate to Galactic globulars, an observational uncertainty of Ac = ±0.2 translates to an 
r.m.s. variation of ±0.2 dex in £ as well. 
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If the mass density profile p(r)/po is replaced by its luminosity equivalent, j(r)/jo, 
the light intensity / at a given projected radius R follows from the standard integral 
1(R) = I{R)/ (j r ) = 2 Jg (j/j ) (r/r ) (r 2 /rg - R 2 /r$y 1/2 dr/r . The central surface density 
Zo is again a function of c only, and it is shown in the bottom left panel of Fig. [l]. The full 
profile T(R) also uniquely predicts the radius TZ(c) = Rh/ro which contains half of a cluster's total 
luminosity in projection. This is plotted against the concentration parameter at the bottom right 
of Fig. Both Xo(c) and 7Z(c) figure prominently in the analysis of §§4 and 5 below. 

It will also be necessary to refer to another feature of single-mass King-model clusters, 
namely, the dimensionless total luminosity C(c) = L/j^r^ = J rt 4tt(j/jo) ( r / r o) 2 dr/rQ. This 
quantity is shown for a range of central concentrations in the top panel of Fig. |2[ The lower half 
then compares the calculated luminosity — -L mo d = £( c )jo r o — to that directly observed for the 
139 globular clusters of Harris (1996). The agreement for the 109 regular clusters is excellent 
(as was also noted by Djorgovski 1993a): (log I/ m od — log L) = —0.017 with an r.m.s. scatter, 
A = 0.25 dex, that is less than the observational errorbar of ±0.3]^ Figure || can therefore be 
viewed as confirming that regular globular clusters are indeed well described by King models — and 
as showing the difficulties with PCC clusters, whose luminosities are often underestimated on the 
basis of their crude King-model parametrizations. 

As has already been suggested, the concentration parameters measured for Galactic globular 
clusters span a reasonably wide range of values, from c ~ 0.5 to c ~ 2.5 (Harris 1996); moreover, 
there is a strong systematic quality to the variation of c from cluster to cluster (see, e.g., Fig. |6] 
below). Thus, Figs. |l] and ^ also serve to illustrate the dynamical and structural non-homology 
that characterizes the Milky Way globular cluster system: The relationships between global and 
local properties of a cluster (i.e., any of the ratios plotted in these Figures) obviously depend 
quantitatively on its internal density and velocity profiles. But the details of these profiles are 
quite sensitive to the King-model parameter c (or, equivalently, Wq), so that even dimensionless 
ratios such as 8 , 1Z, and C can differ significantly between clusters; one globular cannot, in general, 
be turned into another simply by applying a single scaling factor to its basic physical parameters. 

This is a natural consequence of the fact that globular clusters are well described by King 
models (it would not be the case if, for example, these objects were singular isothermal spheres), 
and it is clearly key to a full understanding of cluster structures and a complete explanation of the 
various correlations between their properties. Although this point has always been appreciated 
in principle (e.g., Djorgovski 1995; Bellazzini 1998), its consequences are fully explored for the 
first time in this paper. Throughout, the functions of c just discussed are used repeatedly to 
derive and manipulate physical parameters of the Milky Way globulars, thus always taking into 
account — or, in some sense, correcting for — any non-homology. Whenever an evaluation of £, Iq, 



4 This errorbar, drawn in the bottom panel of Fig. ^, follows from (conservatively) combining in quadrature the 
uncertainties of ±0.1 in log ro; ±0.12 in log (jofo)', and ±0.16 dex in log C(c) (which follows from the shape of C in 
the interval 0.5 < c < 2.5, given Ac = ±0.2 dex). 
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1Z, or C is required, it is obtained from an observed c-value by interpolating on the numerical 
model curves plotted in Figs. |l] and ^. (For convenience, however, Appendix B also gives simple 
approximations to these dimensionless quantities as polynomial functions of the concentration c.) 
It is worth noting that analogous corrections for non-homology in elliptical galaxies may remove 
some significant part of the observed "tilt" of their fundamental plane (see Graham &; Colless 
1997; Busarello et al. 1997). 

3. Two Basic Properties of Galactic Globular Clusters 

The full spatial structure and internal dynamics of a single-mass, isotropic King (1966) model 
cluster are set by specifying four independent physical parameters. The four most convenient, 
for purposes of comparison with observations of real clusters, are the central concentration 
c [or the central potential depth Wo, to fix the overall shape of a surface brightness profile 
u(R) ~ —2.5 log I(R)]; the scale radius ro [for horizontal normalization in a plot of fi(R) 
vs. projected radius R]; the central luminosity jo (the vertical normalization, related to /i as 
described above and in Appendix A); and the core mass-to-light ratio To (for normalization of 
the model to observed line-of-sight velocity dispersions). More generally, however, any set of four 
linearly independent combinations of these quantities will serve equally well as a physical "basis" 
for the sequence of King models. (Such additional factors as metallicity or galactocentric position 
are quite separate from any model characterization of a cluster's internal structure, and they are 
best viewed as "external" physical parameters.) As will be described below, the basis ultimately 
chosen here comprises the concentration, c; the (logarithmic) mass-to- light ratio, log To; total 
luminosity, log L (which agrees well with the model quantity [log jo + 3 log ro + log £(c)]); and 
binding energy, log Ef, = const. + [2 log To + 2 log jo + 5 log ro + log £(c)]. 

At any rate, the implication is that, insofar as Galactic globulars can be approximated by 
the simplest King models — that is, to the extent that internal velocity anisotropies and a range of 
constituent stellar masses can be ignored as second-order corrections — they constitute a nominally 
four-parameter family of stellar systems. As was mentioned in §1, however, there exist many 
correlations between various properties of globular clusters in the Milky Way, implying that these 
objects in reality inhabit only a small part of the physical space potentially available to them in 
principle. Indeed, as Djorgovski (1995) and Bellazzini (1998) have suggested, and as will now be 
shown in a different manner, there exist at least two (and perhaps three) independent constraints 
linking the four basic parameters of globular clusters, so that they are confined to a fundamental 
plane (perhaps a line) in the theoretical 4-space of King models. 
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3.1. A Constant Core Mass-to-Light Ratio 

As was discussed in §2, the central line-of-sight velocity dispersions a Pt o compiled for 39 regular 
and 17 PCC globular clusters by Pryor & Meylan (1993) have been converted to King-model scale 
velocities o"o using the concentration parameters c given by Harris (1996). With measurements of 
all three of <7o, ro, and jo thus in hand, self-consistent mass-to-light ratios Tyo may be calculated 
for fully one-third of the Milky Way globulars via the core-fitting procedure of Richstone &; 
Tremaine (1986): Given do and ro for any one cluster, equation @ yields a core mass density po, 
from which follows T^o = Po/jo- 

Figure ||] plots po against jo for the Pryor & Meylan dataset, with the regular clusters 
represented by circles (open symbols use the directly observed a Pt o in place of the corrected 
do in eq. |J) and PCC objects shown as open squares. The straight line traces the relation 
log po = 0.16 + log jo, obtained by a least-squares fit to the non-core-collapsed clusters (filled 
circles) only. Clearly, the intercept in this equation is the mean mass-to-light ratio (log Tyo). 
Moreover, the r.m.s. scatter about this fit (A = 0.2 dex) can be attributed entirely to random 
observational uncertainties;^] there is no evidence for any significant variation in Tyo from one 
cluster to another. (The apparent departure of the densest core-collapsed clusters from the 
mean line in Fig. || is certainly not entirely real, and it may be completely spurious: forcing a 
King-model fit to these objects leads to underestimates of their core densities jo, and thus to the 
underestimates of their total luminosities in Fig. || and to overestimates of Tyo here. Once again, 
PCC clusters are never included in any of the quantitative analyses in this paper.) 

This result has two implications. First, almost as a practical matter, it can confidently be 
assumed, even in the absence of direct velocity-dispersion measurements, that all other (non-PCC) 
Galactic globulars share the same mass-to-light ratio, viz. 

log Tyo = log (po/jo) = 0.16 ± 0.03 [M Lq 1 ] , (4) 

where the uncertainty is just 0.2/\/39. This mean Tyo = (1.45 db 0.1) MqLq 1 is the same as 
the direct average of the values given for individual clusters by Pryor & Meylan (1993), even 
though these authors define the mass-to-light ratio differently (as the ratio of central surface 
densities) and fit different (multi-mass) King models to the raw velocity and surface-brightness 
data. Their numbers also show no evidence for any significant dependence of Tyo on jo, or on 
cluster luminosity. Figure |3| is also roughly consistent with the mean Tyo — 1.2 MqLq 1 found 
by Mandushev, Spassova, & Staneva (1991) (who, however, claim marginal evidence for a weak 
trend, Tyo oc £ - 08 ± - 06 with large scatter). 

The second consequence of Fig. ^ and equation ([|) is an effective reduction of the 
dimensionality of the parameter space available to Galactic globular clusters, from four to three; 



5 Given A logoo — ±0.09, A log ro = ±0.1, and A log (joro) = 0.12, and with Yy,o oc cro/jo^o the uncertainty 
A log Tv,o is typically ±0.24 dex even if the constituent errors are uncorrelated. 
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in practice, the core mass-to-light ratio is not a free physical variable. This point is taken up 
again in §4.2 below, where it is shown that the constraint log Ty,o = 0.16 is equivalent to one of 
Djorgovski's (1995) bivariate correlations defining the globular cluster fundamental plane. 

3.2. Binding Energy vs. Luminosity 

Following the discussion at the beginning of this Section, the King-model parameters 
remaining to characterize Galactic globulars are c, ro, and jo — or any three linearly independent 
combinations of these quantities plus a constant Ty,o- The remainder of this Section, as well 
as §4, is devoted to establishing total luminosity L (or mass Ty^L) and binding energy E b as 
replacements for ro and jo as two of the axes in the parameter space of King models. There 
are several justifications for this replacement. First, there exists a previously unrecognized, 
but exceptionally strong, correlation between E b and L in the Milky Way cluster system. This 
correlation further reduces the dimensionality of the space inhabited by real globulars, and it 
provides an immediate physical explanation for Djorgovski's (1995) second equation for the 
fundamental plane (§4.2). Second, an emphasis on binding energy in particular facilitates 
comparison with the expected properties of gaseous protoclusters (McLaughlin, in preparation), 
i.e., this quantity is a natural focus for theoretical investigations of the origin and evolution of 
the globular cluster fundamental plane. And third, a view in terms of binding energy might 
profitably be extended to galaxies and clusters of galaxies, perhaps to yield new insight into their 
fundamental planes as well. 

Given equations @ and (||), E b may be calculated for a King-model globular cluster using 
measurements of c and any two of a$ (which, again, is inferred from the directly observed a Pt Q and 
c), ro, and po = Ty.oJo in one of the three equivalent relations, 

E b (a ,r ) = 4.639 x 10 45 erg (^pr)' (^) £(c) (|a) 

E b (a , P0 ) = 5.995 X 10- erg (^) 5 (^) ^ 8(c) <g» 

E b (r , P0 ) = 1.663 x 10 41 erg (^g^f S(c) . #) 

With Ty.o a constant, equation (j3]c) is of particular interest in that it allows for an estimate 
of binding energy based entirely on photometry, i.e., velocity-dispersion measurements are not 
explicitly required. The function S (c) is shown in Fig. [l| above, and it can be approximated by 
equation (|B2| ) below. As was mentioned in §2, an r.m.s. variation of ±0.2 dex in c translates to 
A log 6(c) = ±0.2; thus, the typical measurement errors in c, log ro, and log jo result in a total 
observational uncertainty of ±0.45-0.75 dex (depending on whether or not the individual errors 
are correlated) in the logarithm of E b . 
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Figure ||] plots the concentrations and the binding energies, calculated using each of equations 
(|5|) (supplemented by eq. Ej), against luminosity and mass for the 39 King-model globulars 
(circles) and 17 PCC clusters (squares) with eg known from Pryor &; Meylan (1993). Note that 
the total masses used in this Figure are taken directly from Pryor & Meylan, who derived them 
from fits of anisotropic, multi-mass King models to the cluster data. These masses are therefore 
highly model-dependent, and in fact it is not clear what relevance they have in the analysis of 
this paper, which treats globulars as isotropic, single-mass clusters (in particular, the equality 
M = TyqL is implicit here, but it does not necessarily hold for the numbers of Pryor & Meylan 
1993). Nevertheless, the right column of Fig. |I] demonstrates that correlations between cluster 
concentration and luminosity/mass, and between E^ and L or M, exist independently of specific 
model assumptions. 

The top row of Fig. || shows a rough trend of increasing central concentration with increasing 
cluster luminosity. Such a correlation has also been noted by, e.g., van den Bergh (1994), 
Djorgovski & Meylan (1994), and Bellazzini et al. (1996), and it was recognized in essence long 
ago, by Shapley & Sawyer (1927). This dependence is visibly weaker than that of E^ on L, an 
impression that is quantified by the difference in the Spearman rank correlation coefficients s for 
the various relationships. (These nonparametric measures of correlation strength, specified in each 
panel of Fig. |], are computed for the regular, King-model clusters [filled circles] only.) A linear 
dependence can be fit to these data — c ~ const. + 0.4 log(L/L ) — but the r.m.s. scatter about 
it is A = 0.3, significantly larger than the observational errorbar Ac = ±0.2 dex. Clearly, c is 
not completeley independent of log L, but the relation between the two does not appear to be 
one-to-one (or, if such a strong link does exist, it is nonlinear). 

On the other hand, log E^ and log L are very strongly correlated. Indeed, if the correlation 
coefficients s in the lower six panels of Fig. |I| are compared to those given by Djorgovski & Meylan 
(1994), Djorgovski (1995), and Bellazzini (1998) for a large number of empirical relations between 
other observables of Milky Way clusters, it can quickly be seen that the correlation between binding 
energy and luminosity (or mass) is as strong as or stronger than any other correlation between 
properties of Galactic globulars. In addition to this, it is clear that a simple power-law dependence 
of Ef, on L or M completely describes the data. Least-squares fits (for the regular clusters only) 
of log Ef, against either log L or the log M from Pryor & Meylan (1993) are drawn in Fig. ||. The 
r.m.s. scatter of log Ef, about the straight lines in these six fits always lies within the observational 
errorbar (±0.5 dex) shown in the second row of the Figure. 

The three fits of log Et, against log L are statistically indistinguishable. This is only to 
be expected; it reconfirms the main underlying assumption here — that globular clusters are 
well described by isotropic, single-mass King models — and the finding of §3.1, that the core 
mass-to-light ratio is an essentially constant log Tv,o = 0.16. Figure || specifically compares the 
binding energies computed using measurements of ro and po (eq. [j|c]) to those calculated from gq 
and ro (eq. ||a]). Since the ratio of these is just (inGTyfljorQ /9ctq) 2 , it should be — and is — equal 
to unity in the mean, with an r.m.s. scatter (A = 0.4 dex) just twice that found in Fig. |3| for the 
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basic identity Ty,o = Po/jo- 

If log Ej, is then specified, for each of the 39 regular clusters in Pryor & Meylan's (1993) 
sample, as the direct mean of the three values derived from equations (|5|), a least-squares regression 
against cluster luminosity results in 

log (Eb/exg) = (39.82 ± 0.75) + (2.05 ± 0.15) log (L/L Q ) , (6) 

with a Spearman rank correlation coefficient s = 0.92 and an r.m.s. scatter in log Ef, of ±0.43 dex. 
This relation is significantly different (the uncertainties in eq. || are ±1<t) from the Eb oc L 1 ' 5 
proposed by Saito (1979). 

It is perhaps worth noting that the three fits of log Eb vs. log M in Fig. || are also statistically 
identical to each other, but that they differ formally from the Eb(L) fit: Eb oc M 1 - 8±0 - 1 is indicated, 
as opposed to Ef, oc £ 2 - 05 ± - 15 . At first glance, the difference in these exponents might seem to 
imply a varying global mass-to-light ratio, (M/L)y oc ^ 014±011 Q r so. Such a result is obviously of 
very low statistical significance, however, and of dubious origin besides. It could well be completely 
an artifact of a comparison between different King models, namely, the single-mass and isotropic 
version used here to compute Eb, vs. the multi-mass and anisotropic variant used by Pryor & 
Meylan (1993) to obtain M. In addition, a more direct inspection of the Pryor & Meylan masses 
shows no clear evidence for any systematic variation of global (M/L)y with cluster luminosity 
(consistent with the lack of any such trend in Ty,o from Fig. ^). As a whole, then, the current 
data and models still show no conclusive signs of any variation in global mass-to-light ratios 
among Milky Way globulars. The rest of this paper will therefore discuss cluster binding energies 
explicitly as a function of luminosity only, and it will be taken as given that Eb scales with total 
mass in the same way. 

Figure [6| shows the concentration c and the binding energy Eb (calculated from eq. ||c] with 
log Tyo = 0.16) as functions of luminosity for all of the 139 globulars with data taken from Harris 
(1996). As before, the 109 regular clusters are plotted as filled circles, and the 30 PCC objects 
as open squares. Evidently, the correlations seen in Fig. |4| extend to the whole of the Milky Way 
cluster system. In particular, the power-law nature and high statistical significance of the Eb vs. L 
scaling are unchanged from the smaller sample of Pryor & Meylan. It is now found that 

log (Eb/erg) = (39.89 ± 0.38) + (2.05 ± 0.08) log (L/Lq) , (7) 

where the uncertainties are again ±1<t estimates. The r.m.s. scatter about this fit (to the regular 
clusters only) is A = 0.53 dex in log Eb, consistent with the combination of purely random 
measurement errors in c, log ro, and log jo. With a Spearman rank statistic of s = 0.93, the 
correlation between Eb and L is better by far than those between any other cluster properties in 
this expanded (and now essentially complete) dataset. 

The rough correspondence between c and log L seen in Fig. |] also persists for regular, 
King-model clusters — c ~ const. + 0.4 log L is again indicated — but this remains less significant 
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(smaller s) and more scattered (the r.m.s. deviation deviation from the straight line in the top of 
Fig. | is A ~ 0.35 dex, again exceeding the typical observational uncertainty in c) than the E^-L 
relation. 

The cluster luminosities plotted in Figs. |] and |6] are the integrated absolute magnitudes 
tabulated by Harris (1996), and they are independent of any King- model fits; they have not been 
computed from the formula L = £(c)jorg. Nevertheless, Fig. |2| above demonstrated that the model 
luminosities do correspond very closely to the directly observed values. Thus, a point of potential 
concern here is that a plot of Eb against L is effectively one of (jg^o) a g am st (jo^o), raising the 
possibility that any correlation might be either trivial or spurious. Figure |7| illustrates that the 
situation is rather more subtle than this in reality, and that the Ef,(L) dependence found here is 
indeed nontrivial. The top two panels of this Figure show log rg and log (j'o r o) (which again is 
closely related to the central surface brightness fiv,o) vs - total luminosity for all of Harris' (1996) 
clusters. It is immediately obvious that ro and joro each correlate with L — the dashed lines drawn 
describe ro oc L~ ' 3 and joro oc L 1 ' 25 — but the r.m.s. scatter A in these plots is much larger than 
the observational errorbars (drawn in the upper left corners) in either case. The point of physical 
interest is that the scatter in log ro is (anti)correlated with that in log (jo r o)- One example of this 
has already been shown in Fig. ^, with its plot of £(c)j'orQ against L; another is now seen in the 
bottom panel of Fig. which shows explicitly that the specific combination (Jo r o) ^ Eb/£(c) is 
much more tightly correlated with cluster luminosity [the errorbar in the upper left, ±0.4 dex, 
represents the quadrature sum of the uncertainties in log (jo r o) an d log Tq\. This already suggests 
that the relation between Eb and L is the more fundamental one, and that the (degraded) vq-L 
and jorQ-L correlations are properly viewed as deriving from it. 

The details of this derivation are discussed in §5 and Appendix A, along with other examples 
of similarly weak or scattered monovariate structural correlations in the Galactic globular cluster 
system. A key consideration is the effect of the clusters' structural non-homology (see §2). 
Thus — for example — the bottom of Fig. [7] is still not equivalent to the lower half of Fig. p], because 
j'oro differs from Eb by a non-constant factor of £(c) (eq. []|c]). Indeed, the straight line in Fig. [j] 
represents the scaling j^r^ oc L 1 ' 7 ; it is only because log £ ~ 0.85c (very roughly, from Fig. |l| 
above) and c ~ 0.4 log L that Eb oc L 2,05 obtains in the end. Note also that the r.m.s. scatter A 
in log Eb vs. log L is actually slightly smaller than that in log (jo r o) vs - 1°6 despite the large 
scatter in the c-log L correlation; this is further evidence that the Eb(L) relation is nontrivial and 
physically significant. 

There are other ways to see this dependence, however, and it might have been anticipated on 
the basis of some previous work. In particular, Fig. |8] shows the projected half-light radii of the 
clusters in Harris (1996) as a function of their total luminosities and (3D) Galactocentric radii r gc . 
This comparison, which is not new (cf. van den Bergh et al. 1991; van den Bergh 1995), shows 
that depends only weakly on L but tends to increase systematically with r gc [the dashed line 
in the upper right panel of this Figure has the equation log (R^/pc) = 0.23 + 0.4 log (r gc /kpc)]. 
Moreover, it turns out that the latter variation is responsible for what little dependence Rh does 
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appear to have on L: as Fig. || also shows (see also van den Bergh et al. 1991), the normalized 
quantity R* h = Rh(r gc /8kpc)~ 0A is essentially independent of cluster luminosity.^] In fact, with a 
rank correlation coefficient of s ~ —0.1, R h and L are about as close to being perfectly tmcorrelated 
(s = 0) as Ef, and L are to being perfectly correlated (s = 1). 

The point of this is that E\, oc JqTqS (c), L = jorQjO(c), and Rh = rolZ(c) together give 
E h oc {£K/C 2 )L 2 /R h . But it happens that the non- homology factor (S1Z/C 2 ) varies very little 
over the range of concentrations 0.5 < c < 2.5 appropriate for regular Galactic globulars (see 
Figs, [l] and [2] above, or Fig. [l4] below). The simple scaling E^ oc L 2 / 1 Rh, which is self-evident if the 
question of non-homology is ignored altogether, is then quite accurate; and the weak dependence 
of Rh on L — even if the influence of Galactocentric radius is also ignored — leads inevitably to 
the rough expectation, Ej, oc L 2 . That is, Figs. || and || are equivalent characterizations of the 
Milky Way globular clusters; Ej, ~ L 2 if and only if Rh ~ constant for King-model clusters, and 
conceptually one constraint is no more correct or fundamental than the other. As a practical 
matter, however, note that the r.m.s. scatter about the mean log R^ in Fig. |§] is about twice the 
measurement uncertainty of ±0.1 dex, whereas the scatter about the line E^ oc L 2m in Fig. [6| is 
basically comparable to the observational errorbar. The source of the scatter in Fig. [8| is discussed 

, which also considers another known correlation (o"o oc L 0,52 ; see Fig. ^ and 
Djorgovski & Meylan 1994) that is equivalent to the binding energy-luminosity scaling found here. 

The rest of this paper keeps its focus on binding energy, rather than half-light radius or 
velocity dispersion explicitly, as one of the four defining parameters of globular clusters; and the 
Ei{L) relation, rather than the constancy of Rh, will be used (along with the constancy of Ty,o) 
to define the fundamental plane. Figure ^ makes it clear that this course is not strictly necessary, 
but it has been adopted for the convenience and insight it offers (as has already been suggested, 
a discussion in terms of Eb can greatly facilitate theoretical studies of cluster formation and 
evolution) . 

Before the globular cluster fundamental plane can be constructed explicitly, however, Fig. [8| 
raises another issue to be addressed: Given that £5 oc L 2 /Rh and Rh oc rg C 4 , it is clear that 
"secondary" quantities such as Galactocentric position or cluster metallicity — which are extraneous 
to the formal King-model fitting process — are not completely disconnected from the basic cluster 
attributes Ty,0; c, L, and E^. What, then, is the total effect of such environmental factors on 
these fundamental parameters and the relationships between them? 



6 The relation claimed by Ostriker & Gnedin (1997), Rh oc L~ 0,63 for clusters with 5 < r gc < 40 kpc, is not 
consistent with the earlier analysis of van den Bergh et al. (1991), and it is not confirmed here. Again, the plot of R H 
against L shows that at a fixed Galactocentric radius, cluster sizes are quite insensitive to their total masses. This 
point is important for discussions of cluster destruction timescales, especially in the low-mass regime. 



in §5 (Fig. 17 
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3.3. The Influences of Galactocentric Position and Metallicity 

Figure ^ shows the relative insensitivity of Ty,o, L, and c to r gc and [Fe/H] for globular 
clusters in the Milky Way. The top row incorporates only the sample of Pryor Sz Meylan (1993), 
with Ty^ = 9 CTq/ (4^(7.70 Tq) as in §3.1. Clearly, a cluster's core mass-to-light ratio is quite 
detached from its metallicity. There is perhaps a hint of a slight decrease in Ty.o towards larger 
Galactocentric radii, but this is not statistically significant. Thus, although it would be interesting 
to examine the question more closely by obtaining velocity data for the remaining two-thirds of 
the Galactic globulars, it will be assumed here that Yy,o is completely independent of r gc . (Note 
that the values for the PCC clusters [open squares] are overestimates, due to the underestimation 
of their central densities jo [cf. Fig. 9. Also, the point at r gc ~ 90 kpc corresponds to NGC 2419, 
which appears as an outlier in a plot of o~o against L [see Fig. [l8| below], and for which Ty,o may 
have been underestimated as a result.) 

The middle panels of Fig. ^ illustrate the well known facts (e.g., Djorgovski & Meylan 
1994) that globular cluster luminosities are completely uncorrelated with [Fe/H] — there is no 
mass-metallicity relation — and very weakly anticorrelated with Galactocentric position. However, 
it can be seen that the latter result does not reflect any clear, systematic dependence of L on r gc , 
but appears instead to be driven by a handful of "excess" faint objects at large radii (r gc ^ 10 
kpc) which are not equally represented in the inner parts of the Galaxy. This slight imbalance 
may be due to dynamical evolution from an initial distribution of cluster luminosities which was 
much more insensitive to r gc (see, e.g., McLaughlin & Pudritz 1996), with a larger fraction of the 
low-mass clusters inside the Solar circle having been disrupted by disk shocking and evaporation 
over a Hubble time. Or, alternatively, since some of the outer-halo globulars are somewhat younger 
than average, it could be that clusters formed with lower average masses at lower redshifts (S. van 
den Bergh, private communication). But whatever its cause, this effect does not change the basic 
scalings inferred in Fig. [6| for concentration and binding energy as functions of luminosity. 

Although a plot of c alone against r gc shows a rough anticorrelation (Djorgovski & 
Meylan 1994), this is simply a result of the tendencies for c to decrease towards lower L 
(c rj const. + 0.4 log L, from Figs. || and ||) and for the lowest-luminosity clusters to lie 
preferentially at large Galactocentric radii. The bottom panels of Fig. ^ correct for this, showing 
that at a given cluster luminosity, c is essentially independent of r gc . Equivalently, the mean 
relation between concentration and luminosity does not change with Galactocentric position (see 



also Figs. 12 and below). It is obvious, too, that c is — like Ty,o an d L — impervious to [Fe/H]. 



Figure |l0| shows that similar remarks apply to the scaling, though not to the normalization, of 
the Ef,(L) relation. The top plot does show a tendency for the lowest-energy King- model clusters 
to be found at the largest radii; but this is due in part to the preponderance of faint clusters 
in the outer parts of the Galaxy. If the overall dependence Ef, oc L 2 ' 05 is removed, the influence 



of Galactocentric position alone can be isolated: the straight line in the bottom of Fig. |10| has 



the equation log(^/erg) — 2.05 log(L/L ) = 40.22 — 0.4 log(r gc /kpc). The slope of this line is 
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uncertain by perhaps ±0.1; it is, of course, just that expected from the proportionality Ef, oc I? /Rh 
and the systematic increase of cluster radii with r gc (Fig. ||). Put another way — and as can be 
confirmed either by plotting Ef, r^ 4 against L or by fitting Ef, vs. L directly in a series of narrow 
radial bins — the slope of the Ef,(L) relation in the Milky Way does not vary with Galactocentric 
position, but the normalization decreases systematically towards larger r gc (see also Fig. |l2| ). All 
the data are described completely by 



log (E b /erg) = [(39.86 ± 0.40) - 0.4 log (r gc /8kpc)] + (2.05 ± 0.08) log (L/L Q ) 



(8) 



in which the errobars are ±1<7 and all clusters are normalized to a common r gc = 8 kpc for 
convenience. Note that the r.m.s. scatter about this relation, A = 0.49 dex, is ~10% lower than 
that in equation (|^) and Fig. ||, which ignored the r gc dependence. It can be attributed entirely 
to the effects of random measurement errors in ro, jo, and c, even under the most conservative 
assumption that these uncertainties are uncorrelated. Also, the intercepts in equations (||) and 
(|7|) are consistent with equation (Jsj) given the median Galactocentric radii of the Pryor & Meylan 
(1993) sample of regular clusters (r^ = 8.7 kpc) and the full Harris (1996) dataset (r^ = 7.0 kpc). 



Finally, Fig. 11 shows that metallicity plays no role in setting a cluster's total energy, a result 
which is completely in keeping with the irrelevance of [Fe/H] to any other aspect of the internal 
structure and dynamics of globular clusters (see Djorgovski & Meylan 1994 for further examples). 



3.4. Summary 

To a very good first approximation, Galactic globular clusters can be treated as realizations 
of single-mass, isotropic King (1966) models. A complete basis for the physical description of the 
ensemble is therefore provided by the four quantities log Ty^j c, log L, and log E b . While this 
particular choice of variables is not necessary — any set of four linearly independent King-model 
parameters is permissible, in principle — it is certainly sufficient. Moreover, it has the advantage 
of immediately revealing two strong empirical constraints on the properties of Galactic globulars 
(each of which will, of course, be recovered in any alternative parametrization of King-model 
space): (1) the core mass-to-light ratio is a constant, log (Ty i o/M Lq 1 ) = 0.16 ± 0.03, and (2) 
clusters' binding energies are set by their total luminosities and Galactocentric positions, through 
log Ef, = 39.86 — 0.4 (log r gc /8kpc) + 2.05 log L. As far as the current data can discriminate, 
the Milky Way cluster system appears almost perfectly to obey these two relations — deviations 
from them can be attributed entirely to random measurement errors — and it therefore populates 
a two-dimensional subspace in the four-dimensional volume of King models. The effects of 
metallicity and environment on individual cluster properties are summarized completely (again, 
within the limits of current data) by the r gc -dependence in the normalization of the Ef,{L) relation. 

Equations (Q) and (||) thus define a fundamental plane for Galactic globulars. Various 
representations of this are developed in the next Section. Essentially, central concentration and 
total luminosity or mass are then left to determine the distribution of clusters on the plane; but 
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even this is not random, as there also exists a rough correlation between c and log L (Fig. |6|). 
This dependence is rather weak by comparison with the constraints on log Tyo and log Ef,, and a 
one-to-one c(log L) relation has not yet been found in the data; but, as will also be discussed in §4, 
the situation does suggest the intriguing possibility that Galactic globulars may actually fall along, 
or at least have evolved from, something closer to a "fundamental straight line" (cf. Bellazzini 
1998, who suggested something similar but with a much different physical interpretation). 



4. The Fundamental Plane 

For much of this Section, it will be taken as given that Tyo is a constant for Galactic 
globulars, as in equation (||). This immediately removes the core mass-to- light ratio as a physical 
variable, and the parameter space remaining available to the clusters is only three-dimensional. 

The preceding discussion and Fig. || suggest a physically transparent view of the globular 
cluster fundamental plane (FP) as a thin, tilted slice of (log L, log Ej,, c) space — the volume left by 
the constraint log Tyo = constant — which is seen edge-on along the c-axis (in the log L-log E\, 
plane) and closer to face-on along the Eb axis (in the c-log L plane). As was discussed in 
§3.3, however, the normalization of the E{,(L) relation decreases systematically with increasing 
Galactocentric radius. Thus, a plot of log £7& vs. log L for the entire Milky Way cluster system, 
with its wide range of r gc , actually provides an edge-on view of a collection of separate and distinct 
FPs which intersect the log E\, axis at different points. Figure || can therefore be improved by 
introducing the normalized quantity (cf. eq. ||) 

log El = log E b + 0.4 log (r gc /8 kpc) , (9) 

which simply removes all environmental influences on the FP. 



Figure |12| plots c and log E£ against log L for the 109 regular and 30 PCC clusters catalogued 
by Harris (1996), divided according to whether they lie within the Solar circle (Rq = 8 kpc) or 
outside it. Two points are apparent from this. First, it is confirmed that the r gc dependence in Ef, 
is adequately described by a scaling close to E b oc r^. 4 , and that the mean increase of c with L 
is, although rough, essentially independent of Galactocentric position (cf. §3.3). The solid straight 
lines in both bottom panels follow the relation given by equation (|8|) above; the dashed lines show 
the 3<7 limits on the fitted slope and intercept. Clearly, clusters at r gc < 8 kpc and r gc > 8 kpc 
obey this equation equally well. Similarly, the dashed lines in the top right panel of Fig. [12] are 
linear least-squares fits of c against log L (shallower slope) and of log L against c (steeper slope) 
for the regular clusters at r gc > 8 kpc (excluding Palomar 1, the most obvious outlier). These 
same lines are then re-drawn in the top left panel of the Figure, where they are seen to be equally 
acceptable as crude descriptions of c vs. log L for the King-model clusters with r gc < 8 kpc. 



Second, it can be seen from Fig. 12 that both c and log E% are more strongly correlated with 



log L for clusters beyond the Solar circle than for those within it: the Spearman rank correlation 
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coefficients are higher, and in the case of c vs. log L, the r.m.s. scatter about the linear fits 
is smaller. This point has already been discussed in relation to the central concentrations by 
Bellazzini et al. (1996) and Vesperini (1997), who note that globulars at r gc < 8 kpc have been 
subjected to much stronger dynamical evolution than have those at larger radii (evaporation 
is faster and disk shocks are more severe in denser regions of the Galaxy). Thus, the very 
existence of a c-log L correlation among the more distant clusters suggests that it has a largely 
primordial origin. These authors' quantitative numerical simulations confirm in more detail 
that "evolutionary processes are unlikely to play a dominant role in establishing the c-log M 
correlation" (Vesperini 1997). The same calculations show explicitly that the main effect of the 
evolution within r gc = 8 kpc is to increase the scatter in the final version of an initial c(log L) 
dependence. This is consistent with the appearance of the data in the top half of Fig. |l2|. (It must 
be noted again, however, that the r.m.s. scatter of measured c values about the linear regressions 
at r gc > 8 kpc is still substantially larger than the typical errorbar of ±0.2 dex. It remains unclear 
whether there might exist a nonlinear c-log L correlation that is significantly tighter, or whether 
there is simply a real scatter in cluster concentrations at any given luminosity.) 

Similarly, the very strong correlation of log E^ with log L for the globulars at r gc > 8 kpc 
implies that this dependence, too, may have been set largely at the time of cluster formation. If 
this is correct, then it seems that the main effect of the stronger evolution within the Solar circle 
may have been a greater fractional depopulation of the low-mass tail of the cluster distribution, 
rather than any significant change in a power-law scaling between Ej, and L. (It is emphasized 
again that the r.m.s. scatter in the data of both lower panels of Fig. ^ is essentially the same as, 
or even smaller than, the typical observational errorbar on log Ej,.) Indeed, an initial scaling close 
to Eb oc L 2 can be expected to be roughly preserved during the course of dynamical evolution of 
an entire cluster system. Direct calculations (see the discussion and references in Murray & Lin 
1992; Harris & Pudritz 1994) show that the half-mass radii of individual clusters are fairly well 
preserved during dynamical evolution, even over the course of a Hubble time, i.e., dRh/dt ~ 0; 
but this and Ef, oc L 2 / together imply d log E^/dt ~ 2d log L/dt. Thus, as a cluster loses 
mass to evaporation and tidal shocking it should move down and to the left in the (log L, log 
plane, along a line roughly parallel to the initial correlation (if this was indeed roughly ~ L 2 , and 
assuming that any changes in 0.4 log r gc can be neglected). 

To put this another way, it is plausible that clusters remain more or less in the fundamental 
plane as they evolve and, thus, that globular clusters may have appeared on the FP as a result of 
the formation process. In addition, the relation between c and log L implies that the distribution 
of globulars over the FP is non-random and also reflects something of initial conditions. To see 
this properly, however, requires a further manipulation of the basic cluster variables. 
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4.1. e-Space 

By construction, the log L, log E£, and c axes are mutually orthogonal in the three- 
dimensional cluster parameter space selected by log Ty,o = constant. And in the Milky Way, the 
view of the FP in a plot of log El vs. log L is nearly perfectly edge-on (again, this follows from 
the fact that random measurement errors can account fully for the scatter of observed clusters 
about the relation log E\ = 39.86 + 2.05 log L). The FP is therefore canted at quite a steep angle 
[arctan (2.05) = 64°] to the (log L, c) plane, and the top panels of Figs. |7| and 12 do not offer a 



truly face-on view of it; plots of c against log L suffer from strong projection effects and may not, 
in general, accurately reflect the true dispersion of globulars over the actual surface of the FP. 
(Plots of c vs. log El would be better, but that plane still makes an angle of 26° with the FP and 
projection effects would still be present.) In order to be as rigorous as possible — and hopefully to 
allow eventually for some connection with the fundamental plane (s) of other stellar systems — it 
desirable to correct for this. This is easily done in a procedure analogous to that used by Bender 
et al. (1992; see also Burstein et al. 1997) to define a three-dimensional "/t-space" for elliptical 
galaxies. The present construct for globulars is based specifically on cluster binding energies and 
will therefore be referred to as e-space. 

The easiest way to a face-on view of the FP is to define a new coordinate system by rotating 
the log L— log El plane through 64° (counterclockwise) about the c axis. The new (and still 
mutually orthogonal) axes thus obtained are proportional to 

d = log El - 2.05 log L 

e 2 = 2.05 log El + log L (10) 

e 3 = c 

Alternatively, of course, Galactic globulars are expected, on the basis of equation (g), to lie on 
a line e\ = const.] and the line €2 = const, is perpendicular to this in the original log L-\og El 
plane. In this e-space, then, edge-on views of the globular cluster fundamental plane are obtained 
by plotting ei against either e 2 or £3, while the (€2,63) plane provides the face-on view. 



This is shown in Fig. 13 for the full sample of 139 globular clusters from Harris (1996). The 
two edge-on views of the fundamental plane are drawn as the bold line e\ = 39.86. The scatter 
of the King-model clusters about this line in both the (€1,62) and the (£1,63) planes are just the 
residuals due to measurement errors. The apparent tendency for the 30 PCC clusters to fall below 
the FP is, again, probably spurious: the underestimation of jo leads to low values of E£ and L, 
and thus of ei, for these objects; and £3 = c has been arbitrarily set to 2.5 for many of them. 

The plot of €3 against €2 confirms that Galactic globulars are not distributed uniformly 
throughout their fundamental plane. The correlation of c with log L clearly carries over into one 
between €3 and €2- Moreover, because the clusters do occupy a fairly narrow swath on the FP, 
the potential problems with projection in the c-log L correlation are actually not too severe: 



the dashed line drawn in Fig. 13 has the equation £3 = —(12.5 ± 3.2) + (0.13 ± 0.03) €2, which 
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corresponds (given the definition of £2 and the scaling E£ oc L 2 05 ) to c ~ 0.68 log L + const. 
(cf. Fig. ^). This particular relation follows from fitting only the regular, King-model clusters at 
r gc > 8 kpc, where, again similar to the case for c vs. log L, the correlation between £3 and 62 is 
strongest: The Spearman statistic is s = 0.8 for the globulars outside the Solar circle, and the 
r.m.s. scatter about the best linear fit is A ~ 0.3 dex; but s = 0.4 and A ~ 0.4 dex inside 8 kpc. 

Again taking the view that cluster properties in the outer halo better reflect initial conditions, 
it seems quite clear that Galactic globulars formed in something more like a "fundamental band" 
than a plane per se. Moreover, it is even plausible (but not proven) that they might have been 
born along a fundamental straight line, a one-dimensional locus in e-space. (This is a term coined 
by Bellazzini 1998, who arrived at a similar conclusion from a very different argument. See Figs. [l6| 
and [l9| below for further discussion of Bellazzini's interpretation of the globular cluster FP.) 
That is, it is conceivable that globular clusters formed through a simple process, governed almost 
completely by the mass and Galactocentric position of a gaseous protocluster, which naturally 
resulted not only in Ty,o = const, and E£ oc L , but also in some (unknown) one-to-one relation 
between central concentration and luminosity. 

The viability of this interpretation depends to some extent — similar to what was concluded 
above — on the possibility that there may be some nontrivial function of €3 = c which correlates 
very tightly with £2, i.e., which shows scatter within the observational uncertainties (at least for 
the outer-halo clusters) and which would therefore serve better as the third axis in e-space. On 
the other hand, while globulars at r gc > 8 kpc should have been less influenced by dynamical 
evolution, it is unrealistic to expect that they could have been completely untouched by it. Thus, 
another possibility is that some or all of the "excess" scatter of clusters in the fundamental (e2, £3) 
plane is indeed irreducible, but that it arose during evolution from an initial cluster distribution 
that showed no intrinsic scatter. Theoretical studies of cluster evolution will be key in deciding 
this issue and will be of considerable importance, whatever the outcome, for theories of cluster 
formation. 

It is obviously also of interest to understand the relation between the present e-space and the 
corresponding K-space defined for elliptical galaxies (Bender et al. 1992) and studied for other 
hot stellar systems (Burstein et al. 1997) — or, equivalently, to explain the connection between 
the bivariate correlations for globular clusters (Djorgovski 1995, and §4.2 just below) and the 
analogous relations for ellipticals (Djorgovski & Davis 1987; Dressier et al. 1987; see Pahre, de 
Carvalho, & Djorgovski 1998 for a comprehensive recent discussion). A potentially important 
difference is that one full dimension — the €3 = c axis — of the observational space employed here is 
given over to a parameter describing the broken homology of globulars, whereas the corresponding 
structural factor is incorporated into all three axes of the space used for ellipticals (see Bender 
et al. 1992). While explicit corrections can certainly be made for the effects of non-homology in 
galaxies (Graham Sz Colless 1997; Busarello et al. 1997), these are perhaps less easily pictured in 
K-space than in e-space. The construction of an e-space for galaxies would require the use of a 
single family of structural and dynamical models that is highly accurate in its description of any 
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individual system (in order to estimate total luminosities/masses and global binding energies by 
extrapolation from direct observations of the galaxy cores). 

That real differences exist between the fundamental planes of globular clusters and elliptical 
galaxies is clear from contrasts in certain of their monovariate correlations, or projections of 
their FPs (e.g., Kormendy 1985; Djorgovski 1993b). What remains to be seen is the extent to 
which these disparities might result from fundamentally different Eb(L) relations, as opposed to 
any number of other physical distinctions (examples being a varying core mass-to-light ratio in 
ellipticals [van der Marel 1991]; the presence of dark matter in the galaxies [see Kritsuk 1997]; 
and possibly different systematics in the structural non- homology ) . It may be significant, in this 
regard, that Busarello et al. 1997 find evidence that the kinetic energy per unit mass measured 
within an effective radius in ellipticals scales as (E^/M) oc Oq' 6 or so, whereas the total (E^/M) in 
globular clusters can be roughly — although not perfectly — estimated as {E^/L) oc L L05 ~ Og (see 
Fig. |l^ below for the relation between ctq and L). 



4.2. Bivariate Correlations 

Djorgovski (1995) used a statistical, principal-components analysis to infer the existence of 
two strong, bivariate correlations involving core and half-light properties of Galactic globular 
clusters with measured central velocity dispersions (i.e., those in the catalogue of Pryor & Meylan 
1993): 

log o- p ,o - 0.45 log r = -(0.20 ± 0.01) y, Vja + (4.17 ± 0.2) (0a) 

and 

log <r Pi0 - 0.7 log R h = -(0.24 ± 0.02) (fi V ) h + (4.83 ± 0.26) , (0b) 

where (fj, v )h = 26.362 - 2.5 log (L/2vri^) is the average surface brightness within a projected 
half-light radius; a Pi o is in units of km s _1 ; and ro, Rh are measured in pc. It was the existence 
of these correlations that originally led Djorgovski to conclude that the Galactic globulars are 
confined to a fundamental plane. But the statistical analysis does not, by itself, offer any physical 
insight into the FP. (Djorgovski suggested, as did Bellazzini [1998] after him, that eq. [p4|a] is 
a reflection of the "pure virial theorem" applied to cluster cores. However, as will be discussed 
further below, this is not a fully satisfactory interpretation.) Indeed, note that equations ( jTT| ) 
appear superficially to constitute two constraints on five different observables, and thus it is 
not clear a priori that they need define a plane of any kind. Now, however, with a specific 
physical basis (log Ty,0) c, log Ef,, and log L) having been chosen to describe globular clusters 
in a King-model framework, it is possible to reduce the number of distinct variables in these 
equations to just four, and to show explicitly — by comparing with the bivariate (<ro,ro, fJ>v,o) and 
(o~o,Rh, (fJ'v)h) correlations expected to arise from each of the two empirical relations defining the 
fundamental plane — that equation (|i~I|a) derives from equation (||) above while (|ll]b) is a reflection 
of equation (|6|). Thus, Djorgovski's (1995) correlations together identify the same globular cluster 
FP that has been developed here. 
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The basic definition of the King radius (eq. Q), together with a constant log Ty,o = 0.16 
and the definition of surface brightness, yields a relation between do (not exactly the a Pt o used by 
Djorgovski in establishing his correlations), rg, and /iyo- The details of this derivation are given 
in Appendix A; from equation (A9), 



log (j - 0.5 log r = -0.2 n Vfi - 0.5 log J (c) + 4.241 . (0a) 

Here and throughout, o"o is in units of km s _1 , ro is in pc, and /j,v,o is in mag arcsec -2 . The 
function 2q is the dimensionless central surface density for a King model of concentration c; it is 
given in Fig. |l| above and in equation flB3|) below. An equivalent expression, in terms of cluster 
half-light radii and surface brightnesses, follows from defining (fJ,v)h an d from accounting for 
non- homology in the ratio TZ(c) = Rh/ro (see Fig. |] and eq. |B4| ) ; according to equation ( [All ), 



log a - 0.5 log R h = -0.2 (ji v ) h - 0.5 log [£(c)/K(c)} + 4.640 , f2|b) 
where C is the dimensionless luminosity given in Fig. [2] and equation QB5D - 



Similarly, the analysis in Appendix A (see eqs. [ |A12 | and [ A13| j) shows that the finding 



log (Eb/erg) = 39.82 + 2.05 log(L/L & ) (eq. ||; for clusters in the catalogue of Pryor & Meylan 
1993 specifically) is equivalent to 

log do - 0.775 log r = -0.205 /zy i0 + 0.25 log \c(c) 2m /£ (c) J (c) 2 - 05 ] + 3.943 |L3|a) 

(with the dimensionless energy £ again discussed in §2 and Appendix B), and to 

log cr - 0.775 log R h = -0.205 (fj, v ) h - 0.25 log [£(c)/K(c)] + 4.352 . (|3|b) 

Note that, in general, the r gc dependence in the normalization of the Eb(L) relation (eq. Q) 
should appear as an additional term (—0.1 log r gc + constant) on the right-hand side of equations 
(|l3|). In order to make a more direct comparison with the analysis of Djorgovski (1995), however, 
the effect (which is evidently quite small anyway) has been acknowledged only implicitly here, 
by using a normalization of the E^-L scaling that is appropriate for the median Galactocentric 
position of the globulars in the Pryor & Meylan database. 

Figure |TJ confirms the validity of equations (12 b) and (13b) as a representation of the globular 
cluster FP in terms of half-light quantities. The left panels here show the non- homology functions 
0.5 log (C/1Z) and 0.25 log (£/TZ), as obtained in general from numerical integrations of King 
models (curves) and evaluated specifically for the central concentrations of the 39 regular clusters 
(points) in the sample of Pryor & Meylan (1993). The functions vary only slightly among these 
clusters, and the effects of non-homology happen to be effectively suppressed in this projection of 
the FP. The mean values indicated for each function may therefore be applied to derive nearly 
constant intercepts in equation (|i~2|b) for Tyo = const, and equation (|i~3|b) for E\, oc L 2m . The 



right panels of Fig. 14 then make a direct comparison between the resulting FP "predictions" 



(solid lines) and the Pryor & Meylan data. The agreement is clearly quite good and, as expected, 
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the r.m.s. scatter of the data about the expected correlations is within the range of observational 
uncertainty. 



The top panels of Fig. [14|, and equation fll^b), also represent the basic scaling 
L oc (aft Rh)(C/1Z), which ultimately — given that Ty,o is a constant independent of L — is 
just an expression for the total mass of a King model. Since the ratio C/1Z does not vary 
particularly strongly or systematically with c (or, therefore, with L), this result appears to differ 
significantly from the L oc {a^Rh) ' ' that Schaeffer et al. (1993) find for elliptical galaxies 
and Abell clusters (although note that these authors claim that their scaling holds for Galactic 
globulars as well). This is another indication that the globular cluster fundamental plane may 
indeed differ from those of these other systems. Whether this particular contrast is due to 
fundamentally different trends of Ty;o with L, or to different structural systematics (i.e., to the 
non-homology term C/1Z perhaps being more sensitive to L in the larger systems), is unclear at 
this point. 



The two left panels of Fig. 14 further show the product £1Z/C 2 to be very nearly constant as 
a function of c, and thus they confirm that the expectation Eb oc L 2 / Rh for homologous clusters 
actually is a good approximation for King-model globulars as well (as was discussed in §3.2 above; 
see Fig. ||). 

Figure ^ illustrates the accuracy of equations ([l2]a) and ([u]a) for the globular cluster FP. 
Again, the panels on the left show the appropriate non-homology functions for generic King 
models and for the regular clusters of Pryor & Meylan (1993). The top plot shows (as does Fig. |l|) 
that To is a very weak function of c. Non-homology is therefore a rather minor concern for this FP 
equation as well, and the mean value of Iq can be used to estimate a roughly constant intercept 
for equation (|i"^a); this then compares very favorably to the (non-PCC) data in the top right 
panel of Fig. Il5[ The bottom panels, however, show that things are not always so simple. The 
logarithm of £ 2m /£Iq' 05 in the left panel increases systematically (and close to linearly, as the 
dashed line indicates) with concentration parameter c, and hence with cluster luminosity (and, 
thus, with central surface brightness as well; see eqs. []14|] below). This function is therefore not 
well approximated by a constant, and if it is treated as such — if it is simply assigned its mean 
value in equation (|i~3|a), as was done in the other examples here — the neglected non-homology 
manifests itself as a slope in the observed correlation between (log do — 0.775 log ro) and /iv,o 
which is steeper than the naive expectation. 

Equation (|l3|b) for (log o"o — 0.775 log Rh) vs. {nv)h — a result of Eb oc L 2m — comes closest 
to the basic form of Djorgovski's (1995) equation (|TT|b). The fact that his analysis zeroed in on a 
half-light correlation with the slightly different combination (log a Pt o — 0.7 log Rh) may itself be 
understood as a consequence of the fact that Eb ~ L 2m for the Galactic globulars: since Rh is 
then so nearly constant as a function of luminosity, a correlation between 0.7 log Rh and {fJ>v)h is 
not significantly worse, statistically, than one involving 0.775 log Rh- Meanwhile, equation (Oa) 
for (log o"o — 0.5 log ro) vs. nv,o is clearly equivalent to equation (pT|a) from Djorgovski (1995), 
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which therefore is another statement of the fact that the King-model clusters have an essentially 
constant Ty j0 = (1.45 ± 0.10) M & L~ l . 

Finally, as was noted above, it has been claimed (e.g., Djorgovski 1995; Bellazzini 1998; see 
also Djorgovski & Meylan 1994) that equation (]TT|a) (or eq. [i~2|a]) might follow from globular 
cluster cores, viewed as dynamically distinct entities in their own right, being "virialized" with a 
constant mass-to-light ratio. However, this is not entirely accurate. If cores are defined as those 
parts of clusters within the volume r < ro, direct integrations of King models show that they 
do not satisfy the simplest version of the virial theorem; that is, if kinetic energy is denoted E^, 
then 2Ek(ro) ^ £7&(ro) in general. This is shown in Fig. [l6], which plots the ratio 2Ek/Eb as a 
function of r/ro for three representative King models. [Binding energy is defined, as in eq. (|l]), by 
E b (r) = -(1/2) /J 1 A-kt 2 p4> dr; the kinetic energy is computed from E^{r) = (3/2) 4irr 2 pa 2 dr .] 
This may explain why Bellazzini (1998), who assumed 2Ek/Eb = 1 at r = Tq, inferred a 
mass-to-light ratio of only Ty,o — 0.7 M Q Lq 1 — a factor of two too small — from his analysis of 
Djorgovski's (1995) bivariate correlation for cluster cores. This Section has shown that, while 
Ty,o is indeed a constant in Galactic globulars, the precise form of the correlation between 
(log <7o — 0.5 log ro) and pv,o depends on the definition of ro in equation (||); but although the 
basic scaling there does arise generically from a dimensional analysis of Poisson's equation (or, 
equivalently, from the virial theorem including surface terms which vanish only at the tidal radius), 
the normalization (which helps set the intercept in eq. [11a] or eq. []i~2|a1) is a convenience specific 
to King (1966) models. 



5. Other Correlations 

Sections 3 and 4 have presented the main results of this paper: The binding energies of 
109 regular globular clusters in the Milky Way, calculated within the theoretical context of 
single-mass, isotropic King (1966) models, correlate very tightly with total luminosities, and 
decrease systematically with increasing Galactocentric radius (eq. [^). This result and the fact 
(eq. Q) that the core mass-to-light ratio is a constant (at least for the 39 regular clusters where it 
has been directly measured) then imply the existence of a fundamental plane for globular clusters, 
one which has an immediate and clear interpretation even while accounting for the bivariate 
cluster correlations discovered by Djorgovski (1995). This new physical view of the FP is expected 
to aid in developing theories of cluster formation and evolution. 

Meanwhile, on perhaps a more pragmatic note, it is a mathematical necessity, requiring no 
further proof once a King-model framework has been adopted, that any correlation between any 
set of cluster observables can be obtained simply by treating only r gc , L and (to rather a lesser 
degree) c as independent variables and then manipulating the constraints on £^(L,r gc ) and Tv,o 
according to generic properties of King models. (As was concluded in §3.3, all cluster attributes 
are independent of metallicity, a fact which — although interesting in its own right — justifies the 
neglect of [Fe/H] in this discussion.) No other empirical trend contains any physical information 
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beyond the intrinsic properties of these models and the defining equations @ and (|8|) of the 
globular cluster FP. This point is developed quite generally in Appendix A, and it is now illustrated 
in brief for the Milky Way cluster system specifically. 



Equations (|A14|)-(|AT9|) give expressions for various physical quantities in terms of the basis 



chosen in this paper for King-model parameter space, i.e., as functions of L, c, Ty,0i and Eb, 
allowing for arbitrary Tyo and any power-law scaling E\, = A (L/L ) 7 . Applied to the Milky Way 
cluster system, with log Tyo = 0.16, 7 = 2.05, and log(j4/erg) = 39.86 — 0.4 log(r gc /8kpc), the 
analysis of Appendix A therefore results in the following observable dependences on L, r gc , and c: 

log r = -0.05 log L + 0.4 log (r gc /8kpc) - log (c 2 /£) + 1.681 (0a) 

log R h = -0.05 log L + 0.4 log (r gc /8kpc) - log [c 2 /^ + 1.681 (||b) 

log jo = 1.15 log L - 1.2 log(r gc /8kpc) + log {&/£*) - 5.042 

Myo = -2.75 log L + 2 log (r gc /8kpc) - 2.5 log (£ 3 J /£ 2 ) + 34.766 f|d) 

(/j J y) h = -2.75 log L + 2 log (r gc /8 kpc) - 5 log {c 2 /£v)j + 36.761 (|] 

log a = 0.525 log L - 0.2 log (r gc /8kpc) - 0.5 log {£ / C) - 1.872 ([U 

for cluster radii in pc, luminosities in L@, surface brightnesses in mag arcsec -2 , luminosity 
densities in Lqpc -3 , and velocities in km s _1 . 

Djorgovski & Meylan (1994) present many "monovariate" correlations for Galactic globulars, 
of the type log ro vs. My = 4.83 — 2.5 log L; log do vs. fiyo', log jo vs. log r gc ; and so on. All of 
these can be traced back to some combination of equations (|l4|), and thus to King-model definitions 
and the two fundamental-plane relations. (Additional quantities considered by Djorgovski & 
Meylan include dynamical relaxation times, which are derived from measurements of other 
cluster parameters and thus could easily be included in the analysis here.) Clearly, however, the 
separation of L- and r gc -dependences in many cases might give a skewed or degraded impression 
of the true, underlying physical relationship responsible for a correlation; and so can the neglect 
of structural and dynamical non- homology, i.e., of the increase of (and scatter in) c as a function 
of L, and the concomitant variations in C(c), £(c), TZ(c), and Iq(c) from one cluster to another. 



One example of this is shown in Fig. 17, which follows up the FP prediction of equation (14a 



to improve significantly on the rather poor anticorrelation between scale radius and luminosity that 
was seen in Fig. |7| above for the full set of clusters taken from Harris (1996). It is now clear that 
the only reason any such correlation appears is because the non-homology term log [C(c) 2 /£(c)] 
increases (very roughly) as const. + 0.7 c « const. + 0.3 log L (top panel of Fig. [l?]), far outweighing 
the explicit dependence on log L in equation (|l4|a). The correlation plotted in the bottom panel of 
this Figure takes this and the expected dependence on Galactocentric position fully into account. 
The solid line drawn through the data here is just equation (|i~4|a). The dashed line is a formal 
least-squares fit to the 109 regular-cluster points. Although it has the expected slope of —0.05, its 
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intercept is offset slightly downwards, a discrepancy which traces back to the small — and purely 
observational — bias introduced by assuming that L = C(c)jorQ exactly for the Galactic globulars 
(see the discussion around Fig. [2]). The r.m.s. scatter of the data points about either of these lines 
is a factor of almost 2 smaller than in a plot of log ro against log L alone and appears to be fully 
attributable to random measurement errors.^ 

It is worth noting that this approach also explains the "large" scatter of log values in 
Fig. H above (A = 0.21 dex, as opposed to the observational errorbar of ±0.1 on log Rh alone): 
the scatter about the prediction of equation (|i~4|b) is guaranteed to be the same as that about 



equation (14a), since the former comes from adding log 1Z(c) to each side of the latter. 



All of the other equations above can be similarly dissected, but it is particularly worth 
looking at the correlation of velocity scale with total luminosity (eq. fllffil)- This is compared 
to cluster data from the smaller sample of Pryor & Meylan (1993) in Fig. [l8|. In the top panel, 
non-homology and the expected variation with Galactocentric radius are ignored. A fairly strong 
correlation persists in this case, however, because the ratio £/£ depends only weakly on c (as can 
be inferred from the similar shapes of the individual curves in Figs. [I] and [2] above) and because 
all but one of the regular (non-PCC) Pryor & Meylan clusters lie within 2.6 ^ r gc ^ 29 kpc [so 
that the term 0.2 log (r gc /8kpc) ranges only from —0.1 to +0.1 or so]. The exception to this is 
the bright cluster NGC 2419, marked on Fig. 18, which is located at r gc = 91.5 kpc according to 



Harris (1996). This partly explains its appearance as an outlier relative to the straight line drawn 
through the top panel of the Figure: 

log cr = 0.525 log L - 1.928 , 

which follows from applying an average 0.5 log (£ / C) ~ 0.049 and a median r gc = 8.7 kpc (both 
evaluated using the 39 regular clusters in the sample) to equation (|i~4]f). Clearly, this simple 
version of the full relation already provides a reasonable description of the data (aside from 
NGC 2419), and the monovariate correlation is a rather faithful representation of the FP scaling 
oc L 2 05 and the constancy of Tv,o- 

The situation nevertheless improves somewhat if £ / C and r gc are treated properly as variables. 



This is done in the bottom panel of Fig. |18|, which shows a stronger and more significant (Tq-L 
correlation (i.e., the r.m.s. scatter A is lower and the Spearman rank correlation coefficient s is 
higher), and in which NGC 2419 appears as a less extreme outlier (though an outlier all the same) 
The solid line plotted here traces the full equation (|i~4|f), and it is statistically indistinguishable 
from a formal least-squares fit to the regular clusters (dashed line). The scatter of these data 



7 The observational errorbar on log ro is ±0.1 dex, and an uncertainty of ±0.2 in c translates to about ±0.15 in 
log(£ 2 /£). In addition, the derivation of equation (|l~i|a) in Appendix A involves squaring the equality L — £jo?~o, 
thereby introducing additional scatter of 2 x 0.25/\/T09 ~ 0.05 dex about the mean line (see Fig. ^|). Even added 
conservatively, in quadrature, these three terms lead to an expected r.m.s. scatter of A ~ 0.2 dex, essentially that 
found in the bottom of Fig. h/j|. 
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about either line is again that expected given the typical observational uncertainties]^ 

As was mentioned above, equations ( |i~4| ) can be combined in various ways to produce a large 
number of other correlations. For instance, the sum of ([l4|a) and (|l4|c) gives the full dependence 
underlying the rough correlation between log (jo r o) an d log L seen in Fig. [7| above. Another 
example — and the last to be considered here — is found by eliminating log r gc between equations 
(|l4]a) and (|l4|d) above: 

log r - 0.2 h Vi q = 0.5 log L - 0.5 log (C/2 ) - 5.272 . (15) 



Alternatively, the manipulations leading to equation ( A16 ) below show that this relation is really 
just an expression for the total luminosity of a King model, combined with the definition of surface 
brightness; as such, it is completely independent of any fundmamental plane specifications. 

Equation (pi]) relates to the correlation between scale radius and core surface brightness, 
shown in the top left panel of Fig. [0^, for the 109 regular and 30 PCC clusters from the catalogue 
of Harris (1996). This is of interest because Bellazzini (1998) has claimed that the correlation 
(along with one between ctq and /J-vfl, which can actually be obtained from eqs. [14d] and Jl4|f]) 



results from Galactic globulars having a "constant core mass." However, this conclusion was 
based on the assumption that cores satisfy the virial theorem in the form 2Ek{r$) = Eb(ro), and 
this was shown above (Fig. |i~6| ) not to be true in general. Moreover, the basic hypothesis of a 



constant Mq = Ty,ojo r o ^ s no * borne out by the data. The top right panel of Fig. |19j, for instance 



shows clearly that measured values of j'o r o range over more than two orders of magnitude — far 
in excess of the observational errorbar — in the Milky Way cluster system. Bellazzini (1998) also 
recognizes this but implies, in essence, that the scatter, although unarguably real, is consistent 
with more or less random excursions from a line of constant core mass. However, the bottom left 
panel of Fig. [II] shows that this is not the case either: If it were true that j'o r o ~ const., such that 
r o ~ (jo r o) _1//2 with some intrinsic scatter, then the quantity (log tq — 0.2 fiv,o) would have to be 
uncorrelated with any other cluster property; but instead it shows a significant dependence on the 
total cluster luminosity. 

The point is that both the monovariate ro-//y,o correlation, and a correlation between 
log ro — 0.2fiv,o and log L, are expected on the basis of equation (lilf) . But this says nothing 
about a constant cluster core mass; rather, it implies that King-model clusters at a given total 
luminosity (and with a single central concentration) will always show a correlation of the form 
log ro ~ 0.2/iy,0) by definition. The rougher trends and the large scatter in the top half of Fig. |^ 
reflect the superposition of many such correlations for clusters with a wide range of total L in the 
Milky Way, modified somewhat by systematics in the non-homology terms C and Iq that result 
from the range in c-values as well. 



8 In this case, the errorbar on log oq is ±0.09 dex (§2). 
contribution from 0.5 log Ty.o (see the general eq. |A1S|]), 
Tv.o oc o"o; see §3.1). The total expected scatter in Fig. hq is 



However, the intercept in equation (|l4p) includes a 
which is further subject to errors in log 00 (because 
therefore 2 x A (log a ) = 0.18. 
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The solid line in the bottom left of Fig. 19 is just equation (|T^) with an intercept evaluated 
assuming a constant log (£/!()) = 1-007, the average for Harris' (1996) King-model clusters; the 
dashed line is a least-squares fit to the data. The two differ — and the scatter about either exceeds 
the observational errorbvars on the data — because the ratio C/2q has a significant dependence on 
c, and hence on log L (see Figs. [I] and ||); the cluster non-homology cannot be ignored here. The 
bottom right corner of Fig. [II] therefore shows the full correlation expected on the basis of equation 
(|l~5|) . The agreement with the data is now excellent; the solid, "model" line is indistinguishable 
from a least-squares fit to the regular clusters, and the r.m.s. scatter is fully within the realm of 
random measurement errors. [This last plot is equivalent to one of (/Zjorj]) 1 / 2 against (L/Lq) 1 / 2 , 
so the observed and expected scatter both are just half those in the bottom panel of Fig. ^ above.] 

Thus, the observed correlation between ro and fiyo m globular clusters is controlled entirely 
by the generic behavior of total luminosity in King models, and by the systematic (but scattered) 
increase of c with L found in the Milky Way. There is no need to postulate a constant core 
mass. Indeed, much more generally, in the scheme developed here — where globulars are defined 
completely by L, c, Yyo, and E£ — there is no room for such new constraints not already provided 
by the equations of the fundamental plane. 



6. Summary 

If they are described by single-mass, isotropic King (1966) models, globular clusters are 
fully defined, in general, by specifying just four independent physical parameters. These were 
chosen in this paper to be a mass-to-light ratio (log Yyo), total binding energy (log Eb), central 
concentration [c = log (rt/ro)], and total luminosity (log L). It has been shown that (1) all 39 
regular (non-core-collapsed) Galactic globular clusters with measured core velocity dispersions 
have a common mass-to-light ratio (§3.1), 

log (T Vfi /M Q L-/) = 0.16 ±0.03 ; 

and (2) if it is assumed that this also holds for the other regular clusters in the Milky Way with 
no velocity data, but for which measured King-model structural parameters are available, the full 
cluster system (109 objects in all, excluding core-collapsed members) shows a very well-defined 
dependence of binding energy on total luminosity, modulated by Galacocentric radius (§3.2): 

log (E b /erg) = [(39.86 ± 0.40) - 0.4 log (r gc /8kpc)] + (2.05 ± 0.08) log (L/L & ) . 

The scatter about each of these relations is fully accounted for by the typical observational 
errorbars on log Yyo (— 0.2 dex) and log Ej, (~ 0.5 dex). 

With Yyo essentially a fixed constant and E b known precisely (at a given Galactocentric 
position) as a function of L, only two of the four basic cluster properties are truly independent; 
globular clusters in the Milky Way are confined to a narrow, two-dimensional subregion — a 
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fundamental plane (FP) — in the larger, four-dimensional space of King models. The distribution 
of clusters on the plane is then determined essentially by their c- and L- values. They are not 
scattered randomly, however, as there also exists a correlation between c and log L. Still, this does 
not necessarily mean that the locus of globular clusters is only one-dimensional: a linear c-log L 
relation is of poorer quality than either of the empirical constraints on log Ty,o and log Eb, and 
although there might exist some more complicated function which gives a comparably one-to-one 
dependence of concentration on luminosity, one has yet to be found (§§3 and 4). 

Regardless, the mean trend of c vs. log L is independent of Galactocentric position and cluster 
metallicity, as is the value of log Yy,o and the slope of the correlation between log Ef, and log L 
(§3.3). The normalization of the Ej,(L) scaling is also independent of [Fe/H], but it decreases 
towards larger Galactocentric radii. The equation above accounts for this, and the r gc dependence 
there describes the full extent of environmental influences on the globular cluster FP. 

It was shown in §4 that cluster concentration parameters (and, to a lesser extent, binding 
energies) correlate more tightly with luminosity for globulars at r gc > 8 kpc than for those at 
smaller radii. Thus, since dynamical evolution has likely been less effective at erasing initial 
conditions outside the Solar circle — and since simple arguments indicate that clusters could remain 
on or close to the FP as they evolve — it may well be that most properties of the fundamental 
plane were set more or less at the time of globular cluster formation. Quantitative calculations of 
cluster evolution over a Hubble time in the Galactic potential should be applied to check this. 

A three-dimensional "e-space" was constructed in §4.1 from a simple transformation of the 
cluster parameter space left after specifying log T^o = constant (and after removing the effects of 
Galactocentric radius) in order to obtain a directly face-one view of the FP. It could be interesting 
to attempt a similar construction for elliptical galaxies and clusters of galaxies, i.e., to interpret 
their fundamental planes in terms of binding energy as well, in order to make a more direct 
comparison with the globular cluster FP. 

The equations for Yy,o and for Ej, as a function of L, which define the FP, were shown in 
§4.2 to be equivalent to the two strong bivariate cluster correlations used by Djorgovski (1995) to 
argue for the existence of a fundamental plane for globular clusters. Djorgovski (and Bellazzini 
1998) offered a simple virial-theorem argument (which was shown here to be incomplete) as an 
explanation for one of these correlations, but did not interpret the other. The results of this paper 
have put his results on a firm physical footing. 

Finally, since Tv,o, c, Ef,, and L completely define globular clusters, any correlations 
connecting any of their other properties, or any other trends with Galactocentric radius, 
necessarily derive from generic properties of King (1966) models combined with the empirical 
FP relations and their environmental dependence. Appendix A outlined such derivations quite 
generally, and specialized results for the Milky Way cluster ensemble were presented in §5. 

The picture that has been developed here is a simplification in that it works specifically in the 
context of single-mass and isotropic models for globular clusters. This approximation is evidently 
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an excellent one, but because of it the results of this paper can say nothing about the stellar mass 
functions or possible velocity anisotropy in observed clusters. Nor have the kinematics of clusters 
within the Milky Way (e.g., orbital eccentricities, or bulk rotation of the metal-rich subsystem) 
been considered in any way. These issues aside, however, the preceding discussion may be distilled 
into a small set of four main facts to be explained by theories of cluster formation and evolution: 

• All cluster properties are independent of metallicity. 

• The core mass-to-light ratio, Tv,o = (1.45 ± 0.1) M L , does not vary significantly with 
cluster luminosity or Galactocentric position. 

• The concentration parameter, which controls the shape of a cluster's internal density profile, 
correlates with luminosity. The dependence of c on log L is much more significant outside 
the Solar circle than inside, but even there the relationship is not obviously one-to-one. Its 
basic form is, however, independent of r gc . 

• Binding energy is intrinsically a function of luminosity and is regulated by Galactocentric 
position: Ef, oc L 2m r~® A . Since there is no evidence for significant variations in the global 
mass-to-light ratios of clusters in the Milky Way, -E5 is inferred to scale with total mass in 
the same way as with luminosity. 

Once again, the balance of current evidence suggests that these characteristics of the 
fundamental plane — and, thus, most of the systematics of the Galactic globular cluster system — 
were essentially fixed by the cluster formation process. Moreover, it is not implausible that 
this process was controlled largely by a single intrinsic protocluster parameter — the initial gas 
mass — and adjusted by an external influence depending on Galactocentric radius. Further 
discussion along these lines is left for future work (McLaughlin, in preparation). Regardless of any 
interpretive details, however, it clearly will be important to determine the extent to which the 
globular cluster systems of other galaxies can be described in the same simple terms that apply in 
the Milky Way. 

This work was supported by NASA through grant number HF-1097.01-97A awarded by 
the Space Telescope Science Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., for NASA under contract NAS5-26555. 

A. KING MODELS AND GLOBULAR CLUSTER CORRELATIONS 

The analysis in this paper refers to eleven physical parameters which may be observed or 
derived for globular clusters: the central line-of-sight velocity dispersion <7 Pj o; the scale radius vq, 
which is closely related (but not always identical) to the core radius, at which the luminosity 
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surface density falls to half its central value; the central surface brightness /iyoj which is related to 
the central F-band intensity, Jo, in units of L Q pc~ 2 ; the central luminosity density jo (-L pc~ 3 ); 
the central mass-to-light ratio Tyfl, which gives the central mass density, p = Ty j , and which 
must equal the global mass-to-light ratio in (idealized) clusters of single-mass stars; the projected 
half-light radius Rh and the surface brightness {pv)h averaged within that aperture; the central 
concentration c, which describes the global shape of the cluster's surface brightness profile; the 
total luminosity L; and the cluster binding energy, Eb. 

There are, of course, other quantities of interest (e.g., core and half-mass relaxation times; 
average velocity dispersion within the half-light radius), but within the context of King's (1966) 
model for globular clusters, these other parameters can all be derived from some subset of the 
eleven just listed. Indeed, if the isotropic, single-mass King models are adopted a priori as complete 
descriptions of the internal structures of Galactic globular clusters, then only four of these eleven 
variables are truly independent. These are (1) the mass-to-light ratio Tv,o', (2) the concentration 
parameter c = log (rt/ro), which is directly related to the depth of a cluster's potential well; and 
two other parameters — chosen here to be (3) binding energy and (4) total luminosity — required 
to normalize a dimensionless model to an observed object. The reduction to this physical basis 
requires seven definitions and model relations. (See also Djorgovski 1993a for an outline of this 
procedure, and King 1966 or Binney & Tremaine 1987 for more detailed descriptions of the models 
themselves.) 

First, recall that, where they are available, the directly observed central velocity dispersions 
of Milky Way globulars (Pryor & Meylan 1993) have already been converted to King-model scale 



velocities <7q (in §2 above; see Fig. 1 and eq. |B1|). Also, the so-called core radii tabulated by 



Harris (1996; also Trager et al. 1993 and Djorgovski 1993a) are in fact the ro referred to here as 
scale radii (ro differs from the true projected half-power radius in low-concentration clusters). 

Given these things, the definition of ro (King 1966) yields a cluster's central mass density: 

Tv,oh = 9a 2 /(4irGr 2 ) , (Al) 

where the basic scaling pq oc (oq/GVq) follows simply from dimensional analysis, but the coefficient 
(9 /Air) is specific to the King model. A central luminosity volume density alone is derived from an 
observed ro, c, and central surface brightness, by way of the trivial (model-independent) definition 

Hvja = 26.362 - 2.5 log (J /J© pc~ 2 ) (A2) 

and the model relation 

J = l (c)j r , (A3) 
where Iq is the dimensionless function of concentration shown in the bottom left panel of Fig. [I] and 



approximated by equation flB3|). Equation ( |A1| ) is used in §3.1 to infer the central mass-to-light 
ratio Yy;o for the 39 non-core-collapsed clusters catalogued by Pryor & Meylan (1993). 
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The projected half-light radius, King scale radius, and central concentration are linked 
through the function 

11(c) = R h /r , (A4) 

which is drawn in the bottom right panel of Fig. |l| and approximated in equation ( |B4| ) below. The 
average surface brightness within Rh is then given (for L in Lq and Rh in pc) by 

{Hv) h = 26.362 - 2.5 log (L/2irR 2 h ) , (A5) 

where the total cluster luminosity is given in terms of c, jo and ro by 

L = C(c)j r%. (A6) 



(See Fig. [2] and eq. [B5] for the dimensionless function £.) Finally, the cluster binding energy is 



expressed as (cf. eq. [||) 

tp ^0p n ^V /2 <7q 5 (Ty,oJo)~ 1/2 f^YnfT m2 5cm 

with £ the function in equation ( |B2j ) and the top right panel of Fig. 13. (The last two equalities in 
eq. [A7] follow from the definition |A1|| .) 



To repeat, then, a King-model globular cluster is defined fully by its core mass-to-light 
ratio, central concentration, binding energy, and total luminosity. Real clusters therefore can be 
expected to inhabit a nominally four-dimensional parameter space. However, the dimensionality of 
this space is reduced by one for every independent constraint on the four fundamental quantities. 
Thus, given relations of the type inferred in §3 for the Milky Way clusters, 

Ty,o = constant and Ef, = A (L/L ) 7 , (A8) 

globulars constitute an essentially two-parameter family of objects (although, in principle, Ty,o 
or the Ef)(L) relation may further depend on external, environmental variables; in our Galaxy, 
A oc r~® A according to §3.3 above.) That is, having specified L and c for an object (at some 
given r gc ), its binding energy Ef> and mass-to-light ratio Ty,o follow immediately from equation 
( |A8|) . The product (jo r o) is then given uniquely by equation (A6), and the combination (io r o) 



by equation (A7). This sets the values of jo and ro individually, and all other cluster parameters 



follow from equations ( Al)- flA5| ). If there exists a third basic relation, c = c(log L) — a situation 



which is strongly suggested, although not proven, by current data on the Galactic cluster 
system — then luminosity or mass alone (again, at a given r gc in our Galaxy) controls the globular 
cluster sequence. 

It is clear, as a result of this, that any correlations between any observed or derived cluster 
parameters must trace back to some combination of the equations f|Al|j-(f4^) characterizing King 
models, and/or to relations between the fundamental physical variables T^o, c, Ef,, and L. Several 
examples of such observable correlations are now derived quite generally; their application to the 
Milky Way globular cluster system specifically is discussed in §§4 and 5 above. 



A.l. Bivariate Correlations 



Equation (Al) can be viewed either as a definition of 7*0 or as an expression for T^o i n terms 
of (To, ro, and jo- Either way, it may be re- written as 



(4vrG/9) T Vi oio7*o = ^o 1 5 



or, in appropriate units, 
0"o 



kms 



-l 



2 /„ x -l 

— = 5.9870 x 10~ 3 

PC/ 



T 



M P) L 



0^0 



Joro 



L pc" 



Taking the logarithm of both sides and using equations (|A2| ) and (A3) to introduce the observable 
Hvo then yields 



log (T - 0.5 log r = -0.2/iv,o - 0.5 log J (c) + 0.5 log T Vfi + 4.1610 . 



(A9) 



Given a constant log Ty,o (eq. §) and a very slowly varying 2o(c), this becomes an essentially 
bivariate correlation corresponding to one of Djorgovski's (1995) two equations for the fundamental 
plane of Galactic globular clusters (§4.2). It may also be put in terms of half-light cluster radii 
and average surface brightness by noting a relation between /zyo and {^v)h that follows from 
King-model definitions: Equations (|A2|), (|A3|), (|A6|), (|A4|), and (A5) combined give 



Mv.o 



26.362 - 2.5 log 



lo(c)j r 



L&pc 2 
(fi v )h ~ 2.5 log (2vr) - 2.5 log 



26.362 - 2.5 log 

l (c)K(c) 2 ~ 
£(c) 



2tt 



J (c)^(c) 2 L/2TrR 2 h 
C(c) L Q pc- 2 



(MO) 



and substitution of this in equation (A9) yields 

log (T - 0.5 log R h = -0.2 (fj, v ) h - 0.5 log [£(c)/K(c)] + 0.5 log T Vfi + 4.5601 . (All) 

as another expression of equation (Al). 

A power-law scaling of binding energy with total luminosity leads to another pair of equivalent 
bivariate correlations. With Ef, = A (L/L Q ) 7 as in equation (A8), the first equality in equation 
(|A7D reads 

4r = GA (L/L e ye{ c y l , 

where, again, the normalization A in general may vary with Galactocentric position (or with any 
other environmental factor appropriate to a given dataset). Applying equation ( |A6j ) and inserting 
the usual units, this becomes 
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so that taking the logarithm and using the definition of fxv,o (eqs. |A2| and [jA3|| ) yields 

log a - 0.25 (2 7 - 1) log r = -(0.1 7) Mv,o + 0.25 log [£(c) 7 /£(c) T (c) 7 ] 

+0.25 log ,4 + 2.6362 7 - 11.4166 . (A12) 

And finally, a version that refers to half-light cluster quantities is obtained by applying equations 
flAl|) and flAlOl ): 

log ao - 0.25 (2 7 - 1) log R h = -(0.1 7) (fi v )h ~ 0.25 log [£(c)/K(c)] 

+0.25 log A + 2.83577 - 11.4166 , (A13) 

which is shown in §4.2 to correspond to Djorgovski's (1995) second equation for the globular 
cluster fundamental plane in the Milky Way. 



A. 2. "Monovariate" Correlations 



Given the suite of definitions listed above (eqs. |A1 



), plus the two fundamental-plane 
constraints of equation (A8), any globular cluster observables can be expressed as functions of 
central concentration and total luminosity (and any environmental factors, such as r gc , which in 
general can be subsumed in eq. |A8| ). Some of these relations can appear as monovariate cluster 
correlations if non-homology and variations in c are ignored. 

One such expression follows from the third equality in equation (|A7]) . Together with equations 



(A8) and (A6), this reads 
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or, after taking the logarithm of both sides, 



log r = (2 — 7) log L — log C{cf/S{c) - log A + 2 log T v>0 + 41.2208 



(A14) 



As §5 discusses, this does translate to a rough correlation between r$ and L in the Milky Way (see 
Figs. ^ and |l7|). It also can be used to derive the dependence of Rh = TZ(c) tq on L and c: 



log R h = (2 - 7) log L - log C{cf/8{c)n{c) - log A + 2 log Ty, + 41.2208 , 



(A15) 



accounting for the behavior of observed cluster half-light radii in Fig. || above. (Like ro, Rh is 
measured in pc.) 

and 



Correlations involving surface brightness follow from the basic definitions ( |A2|) , (|A 
(A6), which state that 



log r 



0.2^,0 + 0.5 log (j r 3 ) + 0.5 log J (c) - 5.2724 
0.2 fi Vfi + 0.5 log L - 0.5 log [£(c)/J (c)] - 5.2724 
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for j'o in units of Lq pc -3 and tq in pc. This is itself a bivariate correlation; it is discussed further 
in §5, in connection with a rough correlation between ro and fiyo observed in the Milky Way 
cluster system (see Fig. [H]). But it and equation ( A14| ) also lead to 



log jo = (3 7 - 5) log L + log C(cf/£(c) 



+ 3 log A - 6 log T 



V,0 



123.662 



(A16) 



or to 



fjt V fi = 2.5 (3 - 2 7 ) log L - 2.5 log C(cfl (c)/£(c) 2 - 5 log A + 10 log Ty i0 + 232.466 , (A17) 

which corresponds to the correlation between log jo?*o and log L in Fig. [?] above. In addition, 
equations ( A17 ) and ( A10 ) together, or equations ( |A15| ) and dA5| ) together, imply 



(Hv)h = 2 - 5 ( 3 - 2 7) log L - 5 log C(c) 2 /£(c)K(c) - 5 log A + 10 log Ty )0 + 234.461 . (A18) 



Finally, a correlation between core velocity dispersion and total cluster luminosity can be 
obtained directly from the definition Ei = o" Vo£(c)/G', the empirical scaling Eb/evg = A^L/L^ 1 



and equation (A14): 



log (T = 0.5(7 - 1) log L - 0.5 log [£{c)/C{c)\ + 0.5 log A - 0.5 log Ty, - 21.7218 (A19) 
for c"o in units of km s _1 . This version of the E},(L) relation in the Milky Way is compared to data 



in Fig. 18 above. 



The six quantities on the left-hand sides of equations (|A14 )-( A19 ) are all expressed in 
terms of the King-model basis L, c, Ty,o, and Eb = A(L/L Q )^ (once more, any effects of, say, 
Galactocentric position and metallicity are easily allowed to enter through these variables; see §5). 
Moreover, any correlation between any other set of globular cluster observables may be derived 
from equations ( |A14|) ~( [A19| ) — the results of §A.l, for example, could be thus obtained. As for 
any collection of single- mass, isotropic King- model clusters, there is no new physical content in 
empirical scalings beyond those existing just among L, c, Tyo and Ef,. This follows necessarily 
from the fact that this chosen basis is (by definition) complete. 



B. DIMENSIONLESS FUNCTIONS OF CENTRAL CONCENTRATION 

This Appendix lists ad hoc polynomial fits that describe the variation of certain derived 
quantities in single-mass, isotropic King (1966) models, as functions of the central concentration 
c = log (r^/Vo). These fits have been obtained by comparing with numerical integrations of King 
models with 0.12 <^ c <^ 3.6, or 0.5 < Wq < 16 for Wq the (normalized) central potential. The 
fits generally will not be reliable if extrapolated beyond this range in c, but all Galactic globular 
clusters are described by 0.5 < c < 2.5. 
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The projected, or line-of-sight, velocity dispersion at the center of a model cluster is given 
in terms of the velocity scale parameter ctq — which differs from the true one-dimensional velocity 
dispersion for low-concentration models — by 

log Pi oAx ) = -14.203 log [l + 0.11313 x 10 1 1307 c l , (Bl) 

which deviates from the curve in Fig. [l] above by less than 0.005 dex for c ^ 0.5. 

The dimensionless binding energy (eq. [§] above, and the top right panel of Fig. [l]) may be fit 

by 

log £ = log (GE b /a$ro) = -1.64893 + 2.82056 c + 9.38926 c 2 - 26.0275 c 3 + 30.6474 c 4 

-20.7951 c 5 + 8.61423 c 6 - 2.13978 c 7 + 0.291982 c 8 
-0.0167994 c 9 , (B2) 

with an absolute error of ^ 0.01 dex (relative error <; 1%) for c ^ 0.5. 

The central intensity, or surface density, in terms of the central volume density jo and the 
King radius ro, is (cf. Djorgovski 1993a, and the bottom left panel of Fig. |l| above) 

log J = log (Vioro) = 0.3022 - 7.5726 log [l + 0.2180 x 10 1 1291 c l , (B3) 

to within an absolute error of <^ 0.005 dex in the interval 0.5 <J c <^ 3.6. [Note that in the limit 
c — > oo, this formula gives Iq — > 2.005, slightly lower than the correct value for an isothermal 
sphere, Iq(oo) = 2.018.] For these single-mass models, which have spatially constant mass-to-light 
ratios, Iq is also equal to the dimensionless central mass surface density, £o//°o r o- 

The projected half-light radius, Rh, is related to the King radius, ro = {Qo-q/AkGpo) 1 ! 2 , by 

log K = log (R h /r ) = -0.602395 + 1.36023 c - 1.67086 c 2 + 2.65848 c 3 - 2.71152 c 4 

+1.42555 c 5 - 0.274551 c 6 - 0.0381277 c 7 + 0.0217849 c 8 
-0.00225252 c 9 , (B4) 

which differs by no more than 0.007 dex from the log 1Z obtained by numerically integrating any 
King model with 0.12 ^ c ^ 3.6 (bottom right of Fig. [l] above). 

Finally, the dimensionless total luminosity (top panel of Fig. ^) is given by 

log C = log (L/j rfj = -0.725444 + 1.90743 c + 4.63720 c 2 - 13.0266 c 3 + 14.8724 c 4 

-9.43699 c 5 + 3.39098 c 6 - 0.578169 c 7 - 0.00512814 c 8 
+0.0162534 c 9 -0.00165078 c 10 . (B5) 

The relative error of this expression is 0.5%, in log £, for any 0.5 ^ c 3.6; it may be 
compared to the fit given by Djorgovski (1993a), which applies over a more restricted range in 
central concentration. As with equation (P3|), L also gives the dimensionless mass, Mj po r O' °f a 
single-mass cluster. 
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Fig. 1. — Dimensionless functions of central concentration for single-mass, isotropic King (1966) 
models. Analytic expressions for each are given in Appendix B. The observed concentrations of 
Milky Way globulars lie in the range 0.5 ^ c <^ 2.5. 

Fig. 2. — Top panel: Dimensionless total luminosity C of King models, as a function of central 
concentration. See also equation flB5|), Bottom panel: Comparison of model luminosities £(c)jorjj 
to the directly observed L = 0.4(4.83 — My) for globulars from Harris (1996). Filled circles are 
regular, King-model clusters; open squares are post-core-collapse objects. Straight line is the 
equality log L mo d = log L, and the typical observational errorbar on log L mo( i is shown in the 
upper left corner. 

Fig. 3. — Determination, according to equation (|3|), of the mean core mass-to-light ratio of globular 
clusters in the Milky Way. Filled circles correspond to 39 King-model clusters tabulated by Pryor 
& Meylan (1993); open circles have po computed for these objects using the directly observed a Pt o 
instead of the model velocity scale o~o- Open squares represent PCC clusters. The solid line has the 
equation log po = 0.16 + log jo, obtained from a least-squares fit of the regular-cluster data only. 
The errorbars at upper left show typical observational uncertainties in po and jo- 

Fig. 4. — Concentration parameters and binding energies, as functions of total luminosity, for 
globulars with measured central velocity dispersions (Pryor & Meylan 1993). Filled circles 
correspond to regular, King-model clusters; open squares, to PCC objects. The three lower 
rows present the results of calculating Ef, according to each of equations (||). Integrated absolute 
magnitudes are taken from the catalogue of Harris (1996), while total masses are obtained by Pryor 
& Meylan (1993) from fits of multi-mass, anisotropic King models. 

Fig. 5. — Confirmation that the binding energy computed using equation (|5|c) (with log Ty,o = 0.16) 
is equal to that obtained from (|5]a) for the regular globulars (filled circles) in the catalogue of Pryor 
&; Meylan (1993). Open squares again correspond to the PCC clusters in their sample. This plot 
is equivalent to Fig. [j| 

Fig. 6. — Central concentrations and binding energies of 109 King-model globulars and 30 PCC 
clusters (filled circles and open squares) in the catalogue of Harris (1996). The straight line in 
the top panel traces the rough relation c = —0.57 + 0.4 log L. The solid line in the bottom panel 
represents the least-squares fit log E b = 39.89 + 2.05 log L; dashed lines show the 3a limits (from 
eq. @), log E b = 41.09 + 1.8 log L and log E b = 38.69 + 2.3 log L. 

Fig. 7. — Monovariate correlations log ro vs. log L and log jo r o vs - log L for Galactic globulars in 
the Harris (1996) catalogue, as poor reflections of the more fundamental correlation between log E b 
and log L. Filled circles are regular clusters; open squares are PCC objects. The straight line in 
the top panel has log ro = 1.63 — 0.3 log L; in the middle panel, log joro = —2.86 + 1.25 log L; in 
the bottom panel, log JqTq = —1.33 + 1.7 log L. 

Fig. 8. — Dependence of cluster half-light radius on luminosity and Galactocentric radius. Point 
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types are the same as in Fig. [7|. is seen to increase systematically with r gc (the line in the top 
right panel is log Rh = 0.226 + 0.4 log r gc ), but is essentially independent of L at a given position 
in the Milky Way. R^ is the normalized quantity Rh( r gc/& kpc)~ 04 ; the dashed line in the bottom 
plot has log R* h = 0.59. 

Fig. 9. — Insensitivity of three of the basic King-model cluster parameters to the "environmental" 
factors of Galactocentric position and cluster metallicity. Open squares again represent PCC 
clusters; note that c has been arbitrarily set to 2.5 for many of these. 

Fig. 10. — Dependence of cluster binding energy on Galactocentric radius. Straight line in the 
bottom panel represents the fit E(,L~ 2 ' 05 oc r~® A , as in equation (||), for regular clusters from the 
Harris (1996) catalogue (filled circles). The observational errorbar shown for log is ±0.5 dex. 

Fig. 11. — Independence of binding energy and metallicity for the 139 clusters taken from Harris 
(1996). 

Fig. 12. — Concentrations c and normalized binding energies El = £t(r gc /8kpc)- a4 for regular 
and PCC Galactic globulars within and beyond r gc = 8 kpc. Dashed lines in both top panels 
are c = —1.03 + 0.5 log L and c = —2.44 + 0.8 log L. The r.m.s. scatter in c is A = 0.29 dex 
(excluding Palomar 1) at r gc > 8 kpc, and A = 0.35 for r gc < 8 kpc. Solid line in the bottom 
panels is the fit log E£ = 39.86 + 2.05 log L, and dashed lines there are the 3a limits (from eq. ||) 
log El = 41.06 + 1.8 log L and log El = 38.66 + 2.3 log L. The r.m.s. scatter in log E% is ±0.5 
dex. 

Fig. 13. — Fundamental plane of Galactic globular clusters (109 regular and 30 core-collapsed) in 
the catalogue of Harris (1996). The orthogonal axes in e-space are defined in equation (|Tl"|). The 
(€1,62) and (£1,63) cross-sections of this volume provide edge-on views of the fundamental plane, 
defined by e\ = 39.86 (eq. ||, drawn as bold solid lines). A face-on view of the fundamental plane 
is found in the (£2,63) slice of e-space. Clusters do not fill this plane uniformly but are confined, 
particularly at large galactocentric distances, to a fairly narrow band within it. The dashed line 
is a fit €3 = —12.5 + 0.13 £2 to the globulars (excluding Palomar 1) with r gc > 8 kpc. Most of the 
PCC clusters appear to fall away from the fundamental plane here, but these objects are generally 
not well described by King models. 

Fig. 14. — Left panels: Effects of cluster non-homology on two empirical, bivariate correlations 
involving half-light parameters of Galactic globulars. The functions C(c), lZ(c), and 8(c) are given 
individually in Figs. |l] and and in Appendix B. Points are placed at the observed c-values of the 
39 regular clusters in the catalogue of Pryor & Meylan (1993). Right panels: Observed correlations 
in the Pryor & Meylan cluster sample. Filled circles are their regular clusters, and the 17 open 
squares are their PCC clusters. Open circles plot the uncorrected <7 Pi o for the regular objects, rather 
than the model velocity scales ctq. Solid lines in the top and bottom panels represent equations 



12|b) and (|l3|b); dashed lines are least-squares fits to the data (filled circles). The distant cluster 



NGC 2419 is noted as an outlier (see Fig. [T 
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Fig. 15. — Left panels: Effects of cluster non-homology on two empirical, bivariate correlations 
involving core parameters of Galactic globulars. The functions Io(c), C(c), and £ (c) are given 
individually in Figs. |l| and |2j, and in Appendix B. The dashed line in the bottom panel shows the 
rough scaling 0.25 log (£ 2 ' 05 /£2q' 05 ) const. + 0.17 c. Right panels: Correlations observed for 



the Pryor & Meylan (1993) clusters. Point types have the same meaning as in Fig. 14. Solid line 



in the top panel represents equation (|l^a) and is indistinguishable from a least-squares fit to the 



regular-cluster data. Solid line in the bottom panel is equation (|l3|a) with an intercept of 4.147, 
estimated by assuming a constant (0.25 log (£ 2m /£lfr 05 )) = 0.204 for all 39 King-model clusters. 
The least-squares fit (dashed line) is quite different (log o~o — 0.775 log ro = —0.246 fAv,o + 4.860) 
because of the significant non-homology in this case. 

Fig. 16. — Virial ratio lE^jE^ as a function of radius inside three single-mass, isotropic King (1966) 
models bracketing the range of central concentration observed for Galactic globulars. Cluster cores, 
defined by r < ro, in general do not obey a simplified virial theorem of the form 2E^ = E^. 

Fig. 17. — Source of the monovariate correlation between log ro and log L in Fig. |7]. Top 
panel: The non-homology term log (£ 2 /£) as a function of central concentration. Solid curve 
is obtained from numerical integrations of King models, points are placed at the observed c- 
values of the 109 regular clusters in the catalogue of Harris (1996), and the dashed line shows 
log (C 2 /£) ~ const. + 0.7 c. Bottom panel: The full correlation involving log ro, c, log L, and 
log r gc . Data for 30 PCC clusters from Harris (1996) are also shown (open squares). Solid line is 



equation (|i4|a): log ro + log (£ /£) = 1.681 — 0.05 [log L — 8 log(r gc /8kpc)]; the dashed line, with 



an intercept of 1.63, is the least-squares fit to the regular-cluster data (filled circles). 

Fig. 18. — Top panel: Monovariate correlation between log o~q and log L for clusters from Pryor &: 



Meylan (1993). Point types have the same meaning as in Fig. 14, The observational errorbar on 
log (Jo is ±0.09 dex. The straight line is equation (|i~4]f) after substituting appropriate averages for 
the ratio £(c)/C(c) and r gc (see text). Bottom panel: The full correlation expected on the basis 
of the fundamental-plane constraints on Tv,o and E£. Solid line is equation (fL"4|f); dashed line is a 
least-squares fit to the regular-cluster data (filled circles). 

Fig. 19. — Top left: Monovariate correlation between log ro and fiv,o for globulars in the catalogue 
of Harris (1996). Solid points refer to the 109 King-model clusters; open squares, to the 30 PCC 
objects. The observed r.m.s. scatter in log ro exceeds the observational errorbar of ±0.1 dex. Top 
right: Real variations in the "core mass," Mq oc jo r o- Bottom left: Deviations from the line expected 
if core mass were a constant are not random, but correlated with total luminosity. Bottom right: 
The first three panels are explained as results of the behavior of total luminosity as a function 
of c in King models, combined with the definition of core surface brightness. The straight line is 
equation (|i~5"l). 
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